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ABSTRACT 

Using deep narrow-band (NB) imaging and optical spectroscopy from the Keck telescope and MMT, 
we identify a sample of 20 emission-line galaxies (ELGs) at z = 0.065-0.90 where the weak auroral 
emission line, [O in] A4363, is detected at >3 a. These detections allow us to determine the gas- 
phase metallicity using the “direct” method. With electron temperature measurements, and dust 
attenuation corrections from Balmer decrements, we find that 4 of these low-mass galaxies are ex- 
tremely metal-poor with 12 + log(0/H) < 7.65 or one-tenth solar. Our most metal-deficient galaxy 
has 12 + log(0/H)= 7.24+g;gg (95% confidence), similar to some of the lowest metallicity galaxies 
identified in the local universe. We find that our galaxies are all undergoing significant star formation 
with average specific star formation rate (SFR) of (100 Myr) -1 , and that they have high central SFR 
surface densities (average of 0.5 M 0 yr _1 kpc -2 ). In addition, more than two-thirds of our galaxies 
have between one and four nearby companions within a projected radius of 100 kpc, which we find is 
an excess among star-forming galaxies at z =0.4-0.85. We also find that the gas-phase metallicities 
for a given stellar mass and SFR lie systematically lower than the local M*-Z-(SFR) relation by «0.2 
dex (2cr significance). These results are partly due to selection effects, since galaxies with strong star 
formation and low metallicity are more likely to yield [O III] A4363 detections. Finally, the observed 
higher ionization parameter and high electron density suggest that they are lower redshift analogs to 
typical z > 1 galaxies. 

Subject headings: galaxies: abundances — galaxies: distances and redshifts — galaxies: evolution - 
galaxies: ISM — galaxies: photometry — galaxies: starburst 


1. INTRODUCTION 

The chemical enrichment of galaxies, driven by star for- 
mation and regulated by gas outflows from supernovae 
and inflows from cosmic accretion, is a key process in 
galaxy formation that remains to be fully understood. 
The greatest difficulty in measuring chemical evolution 
across all galaxy populations is the need for rest-frame 
optical spectroscopy. Metallicity determinations can be 
obtained through (1) interstellar absorption lines (e.g., 
Fc II, Mg II), and (2) nebular emission lines (e.g., [O n], 
[O in], and [N Ii] ) . While studies have used absorption 
lines to measure heavy-element abundances (Savaglio et 
al. 2004), the need for deep spectroscopy and complica- 
tions with curve-of-growth analysis have made it diffi- 
cult. As such, the primary method used to measure the 
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metal abundances in galaxies has been nebular emission 
lines. This technique has the advantage of being able to 
probe low-luminosity galaxies since it does not require 
continuum detection. In addition, these emission lines 
can be observed in the optical and near-infrared (near- 
IR) at redshifts of ~3 and below with current ground (see 
e.g., Hayashi et al. 2009; Moustakas et al. 2011; Henry 
et al. 2013; Momcheva et al. 2013) and space-based ca- 
pabilities (see e.g., Atek et al. 2010; Xia et al. 2012), and 
the forthcoming IR capabilities of the James Webb Space 
Telescope (JWST) will extend this further to 2 « 6. 

The most reliable metallicity determination is made 
possible by measuring the flux ratio of the [O ill] A4363 
auroral line against a lower excitation line, such as 
[O ill] A5007. The technique is often called the T e or 
“direct” method for its ability to determine the elec- 
tron temperature ( T e ) of the ionized gas, and hence the 
gas-phase metallicity (see e.g., Aller 1984). However, 
the detection of [O ill] A4363 is difficult, as it is very 
weak (and almost undetectable in metal-rich galaxies). 
For example, the first data release (DR1) of the Sloan 
Digital Sky Survey (SDSS; York et al. 2000) only re- 
vealed 8 new extremely metal-poor galaxies (XMPGs; 
12 + log(0/H) < 7.65) among 250,000 galaxies (Kniazev 
et al. 2003). Even with improved selection and a larger 
sample («530,000), Izotov et al. (2006a) only detected 
[O ill] A4363 at >2er significance in 310 galaxies (i.e., one 
in 1700). While SDSS spectra can be stacked for aver- 
age measurements of [O ill] A4363 (Andrews & Martini 
2013, hereafter AM13), this sacrifices knowledge of the 
intrinsic scatter in the mass-metallicity (Afy-.Z’) relation 
(MZR) , which can also constrain galaxy evolution models 
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(see e.g., Dave et al. 2011). 

While it is unfortunate that the T e method cannot be 
used for the full dynamic range of metallicity, the detec- 
tion in a galaxy of [O III] A4363 alone is a strong indi- 
cation that it is extremely metal deficient. These rare 
XMPGs are suspected to be primeval galaxies that are 
undergoing rapid assembly at the observed redshift (Kni- 
azev et al. 2003, and references therein). Studying larger 
samples of them can provide a better understanding of 
the early stages of galaxy assembly. 

One possibility is that XMPGs have significant out- 
flows that are induced by supernova from massive star 
formation. These outflows could drive large amounts 
of metal-rich gas out of the galaxy, thus decreasing the 
metal abundances. Studies have found that outflows are 
prevalent in metal-poor galaxies through (1) detection 
of outflowing ionized gas from integral field unit (IFU) 
spectroscopy of an XMPG (Izotov et al. 2006b); (2) blue- 
shifted absorption lines from slit spectroscopy of star- 
forming galaxies, suggesting that outflows are ubiquitous 
(e.g., Weiner et al. 2009; Martin et al. 2012); and (3) ev- 
idence that galaxies with higher specific SFR (SFR per 
unit stellar mass; sSFR = SFR/M*) are more metal poor 
(Ellison et al. 2008; Lara-Lopez et al. 2010; Mannucci et 
al. 2010). 

Another explanation is that metal-deficient gas is sup- 
plied either from the circumgalactic medium (CGM) per- 
haps through a “cold-mode” accretion phase (Keres et 
al. 2005; Dekel et al. 2009) or from the strong interac- 
tion with a nearby merging companion (see e.g., Kewley 
et al. 2006; Rupke et al. 2010). 

The two most metal-deficient galaxies known to date 
are I Zw 18 (Searle & Sargent 1972) and SBS0335-052 
(hereafter, SBS0335; Izotov et al. 1990) with oxygen 
abundances of 12 + log(0/H) = 7.14 (Zq/ 35; Izotov et 
al. 2006a) and 7.19-7.34 (Z 0 /( 23-32); Izotov et al. 1999), 
respectively. Efforts have been made to increase the 
galaxy sample with direct metallicity determinations in 
the local universe (Brown et al. 2008; Berg et al. 2012; 
Izotov et al. 2012a), and at higher redshift (Hoyos et 
al. 2005; Kakazu et al. 2007; Hu et al. 2009; Atek et al. 
2011). These studies have either targeted galaxies with 
low luminosity or high equivalent width (EW) emission 
lines. The latter are found using NB imaging, grism spec- 
troscopy, or unusual broad-band colors to select them. 
For example, Kakazu et al. (2007) and Hu et al. (2009) 
utilized NB imaging to select ELGs at z ~ 0.40-0.85 and 
conducted optical follow-up spectroscopy with Keck to 
detect [O ill] A4363 and determine direct metallicity for 
28 galaxies. They found one galaxy where their measured 
metallicity is 12 + log(0/H) = 6.97T0.17. To date, only 
~70 galaxies are known to have 12 + log(0/H) < 7.65 
with the majority (90%) of them at z < 0.1. 

In this paper, we focus on our spectroscopic detections 
of [O hi] A4363 in 20 galaxies at z = 0.065-0.90 (aver- 
age of 2 = 0.54 ± 0.22) in the Subaru Deep Field (SDF; 
Kashikawa et al. 2004). These galaxies were initially se- 
lected for their excess flux in NB and/or intermediate- 
band filters produced by emission lines. In particular, we 
have rest-frame spectral coverage of at least 3700-5010A, 
enabling metallicity determinations using the T e method. 

The outline of the paper is as follows. In Section 2 we 
describe the imaging survey for the SDF, the selection 
of over 9,000 ELGs, the follow-up optical spectroscopy 


that we conducted, and our accurate flux calibration ap- 
proach. We then discuss in Section 3 our approach for 
detecting and measuring nebular emission lines, which 
yields a spectroscopic sample with [O III] A4363 detec- 
tions that are significant at >3o\ We also present ar- 
guments for why all but one of our galaxies is ionized 
primarily by young stars. Section 4 then describes how 
we determine: (1) the dust attenuation properties in 

our galaxies; (2) the electron temperature and the gas- 
phase oxygen metallicity; (3) the dust-corrected SFRs; 
(4) the stellar properties from spectral energy distribu- 
tion (SED) fitting; (5) the nearby environment; and (6) 
the SFR surface density. In Section 5, we compare our 
results to other galaxies that have direct metallicity de- 
terminations, and discuss selection effects for our sam- 
ple. Finally, we estimate the space densities of compact, 
extreme star-forming, metal-poor galaxies found in this 
survey, and consider their implications for the broader 
context of galaxy evolution. 

Throughout this paper, we adopt a flat cosmology with 
Da = 0.7, Dm = 0.3, and H 0 = 70 km s -1 Mpc -1 to de- 
termine distance-dependent measurements, and magni- 
tudes are reported on the AB system (Oke 1974). For ref- 
erence, we adopt 12 + log(O/H) 0 = 8.69 (Allende Prieto 
et al. 2001) for metallicity measurements quoted against 
the solar value, Z Q . Unless otherwise indicated, we 
report 95% confidence measurement uncertainties, and 
“[O ill]” alone refers to the strong 5007A emission line. 

2. THE SUBARU DEEP FIELD 

The SDF has the most sensitive optical imaging in sev- 
eral NB and intermediate (IA) filters in the sky, and is 
further complemented with ultra-deep multi-band imag- 
ing between 1500A and 4.5/rrn. A summary of the ancil- 
lary imaging is available in Ly et al. (2011a) and later in 
Section 4.4. The results for this paper are based on data 
obtained in the NB704, NB711, NB816, NB921, NB973, 
IA598, and IA679 filters. A summary of their proper- 
ties (i.e. , central wavelength, sensitivity) is reported in 
Table 1, and Figure 1 shows the total system response 
through these filters and the surveyed redshifts. 

These SDF data were acquired with Suprime-Cam 
(Miyazaki et al. 2002), the optical imager mounted at 
the prime focus of the Subaru telescope, between 2001 
March and 2007 May. The acquisition and reduction of 
these data have been discussed extensively in Kashikawa 
et al. (2004), Kashikawa et al. (2006), Ly et al. (2007) 
(hereafter L07), and Ly et al. (2012b) for the NB data, 
and in Nagao et al. (2008) for the IA data. In brief, data 
were obtained mostly in photometric conditions with av- 
erage seeing of O'/O-l'.'O for all five NB and two IA fil- 
ters. These data were reduced following standard reduc- 
tion procedures using SDFRED (Yagi et al. 2002; Ouchi 
et al. 2004), a software package designed especially for 
Suprime-Cam data. 

The most prominent emission lines entering these NB 
and IA filters are Ha, [O in], H/3, and [O II], at well- 
defined redshift windows between 2 = 0.01 and 2 = 1.62. 
This results in probing 64% in redshift space and 67% of 
the available comoving volume at z < 1.03. Compared 
to the previous NB survey for XMPGs at 2 = 0.24-0.85 
(Hu et al. 2009), our survey probes 4.7 (3.8) times more 
redshift (volume) space, and is deeper by «1.5 mag in 
the NB filters. 
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Fig. 1. — Top: Total system throughput for the VRqT z' filters (gray dashed lines) and the IA598, IA679, NB704, NB711, NB816, 
NB921, and NB973 filters (black solid lines). Bottom: Redshift ranges probed when Ha, [O III], H /3, and [O II] are redshifted into these 
filters. The SDF probes 64% of redshift space and 67% of comoving volume at z < 1.03. 
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TABLE 1 

Summary of Filters, Emission-line Samples, and Redshift Windows 


Filter 

(1) 

Ac 

[A] 

(2) 

FWHM 

[A] 

(3) 

«11im(3o') 

[mag] 

(4) 

Area 

[arcmin 2 ] 

(5) 

Atotal 

(6) 

Atarget 

(7) 

A spec 
(8) 

z(Ha) 

(9) 

-([O hi]) 
(10) 

*«o a]) 
(11) 

IA598 

6007 

303 

26.79 

870.4 

641 

31 

21 


0.170-0.230 

0.571-0.652 

IA679 

6780 

340 

27.39 

870.4 

790 

76 

54 

0.007-0.059 

0.320-0.388 

0.774-0.865 

NB704 

7046 

100 

26.71 

870.4 

1695 

173 

140 

0.066-0.081 

0.397-0.417 

0.877-0.904 

NB711 

7111 

72 

26.07 

870.4 

1480 

111 

92 

0.078-0.089 

0.413-0.427 

0.898-0.918 

NB816 

8150 

120 

26.90 

870.4 

1602 

300 

204 

0.233-0.251 

0.616-0.640 

1.171-1.203 

NB921 

9196 

132 

26.71 

870.4 

2361 

251 

185 

0.391-0.411 

0.824-0.850 

1.450-1.485 

NB973 

9755 

200 

25.69 

788.7 

1243 

71 

63 

0.471-0.502 

0.928-0.968 

1.591-1.644 

Total 





9264 

870 

713 





Note. — We report the numbers of excess emitters (IVtotal) in Col. (6), the numbers of galaxies with targeted MMT and/or 
Keck spectra (iVtarget) in Col. (7), and the sub-sample with robust spectroscopic redshift (N 3pec ) in Col. (8). 


2.1. Selection of Emission-line Galaxies 

To select NB and IA excess emitters due to the pres- 
ence of nebular emission line(s), we use the standard 
color excess selection, where photometric fluxes from 
these filters are compared against those for adjacent 
broad-band filters that sample the continuum. Since the 
technique has been extensively used, we briefly summa- 
rize it below, referring readers to Fujita et al. (2003) and 
Ly et al. (2011b). We summarize our excess emitter sam- 
ple in Table 1. 

The measured continuum adjacent to the line is deter- 
mined by the two broad-band filters closest to the narrow 
bandpass. For NB921 and NB973, we start with the z' - 
band. Since the central wavelengths of these filters are 
redder than what the z' filter measures, we correct for 
such differences using the i'-z' color (Ly et al. 2012a, b). 
The remaining filters use a flux-weighted combination 
of either the V- and I?c-band (IA598), the Rc~ and i'- 
band (NB704, NB711, and IA679), or the i' - and z'-band 
(NB816): 

/cont = e/blue T (1 e)/recL (1) 

where /blue and / re d are the flux density in erg s -1 cnr 2 
Hz^ 1 for the bluer and redder broad-band filters, respec- 


tively (e.g., i' and z' for NB816), and e = 0.45 (IA598), 
0.5 (NB704, NB711), 0.75 (IA679), and 0.6 (NB816). 

Photometric measurements for broad-band data are 
obtained by running SExtractor (Bertin & Arnouts 1996) 
in “dual-image” mode, where the respective NB or IA 
image is used as the “detection” image. This works well 
because all broad-band, IA, and NB mosaicked images 
have very similar seeing, so that excess colors are deter- 
mined within the same physical scale of the galaxies. For 
the extraction of fluxes and selection of sources, we use 
a 2" diameter circular aperture. For comparison, the 3<r 
sensitivities for the P, Rc, i ' , and z' data are between 
26.27 and 27.53 mag, which are generally deeper than 
the NB and IA imaging. 

We also exclude sources which fall in regions affected 
by poor coverage and contamination by bright fore- 
ground stars. The unmasked regions cover 870.4 arcmin 2 
for all filters with the exception of NB973, which cov- 
ers 788.7 arcmin 2 . The latter is smaller due to higher 
systematic noise in one of the ten CCDs, which we 
mask to avoid significant spurious detections. In to- 
tal, 123123, 97632, 133273, 119541, 84786, 118097, and 
139585 sources are detected in the unmasked regions of 
the NB704, NB711, NB816, NB921, NB973, IA598, and 
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IA679 mosaics, respectively. Among these sources, 1695, 
1480, 1602, 2361, 1243, 641, and 790 are identified as NB 
or I A excess emitters. 

In certain circumstances, our sources are selected by 
more than one filter. This is due to some fortuitous 
redshift overlap of our NB/IA filters such that differ- 
ent emission lines (e.g., Ha and [O ill]) are detectable at 
the same redshift (see Figure 1). Accounting for dupli- 
cate galaxies, the complete SDF ELG sample consists of 
9264 galaxies mostly at z = 0.01 to z = 1.62 with some 
at higher redshift due to Lya emission. We note that 
the sample presented in this paper supersedes our earlier 
multi-NB studies. Previously, we had selected ELGs in 
four NB filters and examined the evolution of the Ha, 
[O III], and [O il] luminosity functions (L07). A fifth NB 
filter (NB973) was later included, and that filter allowed 
us to study dust attenuation at z = 0.4-0. 5 (Ly et al. 
2012a) and stellar properties out to z = 1.5-1. 6 (Ly et 
al. 2012b). 

2.2. Optical Spectroscopy 

The primary results of this paper are based on op- 
tical spectroscopy taken over the past several years, 
with Keck’s Deep Imaging Multi-Object Spectrograph 
(DEIMOS; Faber et al. 2003) and MMT's Hectospec 
(Fabricant et al. 2005). In total, we obtained 945 optical 
spectra for 870 ELGs, and successfully detected emission 
lines to determine redshift for 713 galaxies or 82% of the 
targeted sample. These spectra were initially obtained 
to confirm that the NB technique efficiently identified 
ELGs. A summary of the spectroscopically confirmed 
excess emitters is provided in Table 1. 

The majority (61%) of our spectroscopic sample was 
obtained from Hectospec. These spectra were obtained 
between 2008 Mar 13 and 2008 Apr 14, utilizing the 270 
mm -1 grating blazed at 5200 A to yield spectral coverage 
of 3650A-9200A. The combination of this grating with a 
fiber diameter of 1"5 yields a spectral resolution of f«6A. 

The typical seeing for these observations was between 
0'.'7 and 1"4, and the data were obtained at airmasses 
below 1.35, with 95% (77%) of the data below 1.3 (1.2). 
Four different fiber configurations were used, with on- 
source integrations varying between 4 and 6 20-min ex- 
posures, which is sufficient for cosmic ray rejection. The 
MMT spectra were reduced following standard proce- 
dures with the iraf Hectospec Reduction Software 12 . 
Since no order-blocking filter was used, the spectra at 
wavelengths longward of 8200A have a significant amount 
of contaminating second-order light (up to 30% for the 
relatively blue flux standard stars). The flux calibration 
in the far-red is therefore unreliable, and we do not use 
it. 

The Keck/DEIMOS observations were obtained in 
2004, 2008, and 2009. The 2004 spectroscopic observa- 
tions have been discussed in Kashikawa et al. (2006) and 
L07, and the more recent data have been discussed in 
Kashikawa et al. (2011). In brief, we constructed 13 slit- 
masks with 1" slit widths. This typically corresponds to 
a spectral resolution of R ~ 3600 at 8500A, and spectral 
dispersion of 0.47A pixel -1 . The 830 line mm -1 grat- 
ing and GG495 order-cut filter were used in all DEIMOS 
observations. This set-up resulted in a spectral cover- 

12 http : //tdc-www .harvard . edu/ instruments/hectospec/. 


age of «5000A-l^m; however the coverage varied along 
the dispersion axis of the slit mask. The typical seeing 
for these observations was 0 , /5-l"0 in 2004 and 0" 7- 1'/ 1 
for 2008-2009 with integration times of 2-3 hr. Almost 
all (86%) of the DEIMOS spectra were obtained at low 
airmasses (< 1 .3) , with the remaining data taken at an 
airmass of 1.35. 

2.2.1. Flux Calibration 

The metallicities that we will determine require mea- 
surements between 3700A and 5010A in the rest-frame, 
thus accurate flux calibration is critical. We follow a 
rigorous approach: we (1) observe spectro-photometric 
standards to account for the wavelength-dependent sen- 
sitivity of each spectrograph; (2) correct for slit losses by 
comparing spectra against broad-band photometric data; 
(3) compare emission-line fluxes against NB photometry 
to assess the accuracy of flux calibration; and (4) com- 
pare spectra obtained on different nights for a few dozen 
multiply-observed galaxies. A more detailed descrip- 
tion of the flux calibration for our MMT /Hectospec and 
Keck/DEIMOS spectra is deferred to Appendix A. In 
brief, our various independent tests and analyses yielded 
consistent results, and demonstrated that the absolute 
flux calibration is reliable at the 0.15-0.17 dex (0.12-0.17 
dex) level for Hectospec (DEIMOS). 

3. THE [O hi] A4363 SAMPLE 

To extract fluxes for strong and weak emission lines 
in these spectra, we fit each line with a Gaussian profile 
using the IDL routine MPFIT (Markwardt 2009). The 
expected location of emission lines was based on a pri- 
ori redshift determined by either the [O ill] or Ha (for 
lower redshift). A local median, (/), is computed within 
a 200A-wide region, excluding regions affected by OH 
skylines and nebular emissions. In addition, the stan- 
dard deviation <r(/) is measured locally. Examples of the 
computed medians and standard deviations are shown 
in Figure 2. To determine the significance of emission 
lines, we integrate the spectrum between l G — 2.5 oq and 
lc + 2.5 <jg, where a G is the Gaussian width: 

+2.5 <tg 

Flux = £ [/(A-+)-(/>]xr. (2) 

— 2.5fTG 

Here, l' is the spectral dispersion (1.2lA pixel -1 for 
MMT and «0.47A pixel -1 for Keck). We then com- 
pute the signal-to-noise (S/N) of the line by dividing the 
integrated flux by: 

Noise = cr(/) x l 1 x + N pixeh (3) 

where Np ixel = 5a G /l'. 

Adopting a minimum significance threshold of 3cr, we 
identify 20 and 14 [O III] A4363 detections with MMT 
and Keck, respectively. We visually inspected each 
[O ill] A4363 detection. For MMT, we found that OH 
sky-lines contaminated [O ill] A4363 in three cases, and 
H f3 in three other galaxies. This results in a final sam- 
ple of 14 [O ill] A4363 detections. For Keck, OH sky- 
lines contaminated [O ill] A4363 in three cases and [O n] 
measurements in two other cases, while three sources 
lack full spectral coverage (missing [O n], H/3, and/or 


Extremely Metal-Poor Emission-line Galaxies 


5 


[O hi] AA4959,5007). 13 This reduced the [O hi] A4363 
Keck sample to 6 galaxies. 

Our final sample of [O ill] A4363 detections, hereafter 
the “[O ill]- A” sample, consists of 20 galaxies. The 
MMT and Keck [O ill] A4363 detections are shown in 
Figure 2 with the full MMT (Keck) spectra provided 
in Figures B1-B2 (Figure B3). We also summarize our 
[O ill] A4363 sample in Table 2. For convenience, our 
galaxies are identified as “MMT” and “Keck” followed 
by a sequential number. In Table 3, we provide the [O il], 
H /3, and [O III] emission- line fluxes along with the O 32 
and R 23 (Pagel et al. 1979) flux ratios: 


[0 hi] AA4959, 5007 
° 32 = [O 11 ] AA3726, 3729 ’ 

D _ [O 11 ] AA3726, 3729 + [O ill] AA4959, 5007 

R23 = TT7, 


(4) 
• (5) 


Also, emission-line luminosities and rest-frame EWs 
(EWo) are illustrated in Figure 3. The latter are de- 
termined by measuring the continuum from the SEDs 
with corrections for emission-line contamination (see Sec- 
tion 4.4). The former is L = 47rd 2 Thine, where <#, is the 
luminosity distance, and TLine is the emission-line flux. 

We note that for two galaxies (MMT#04 and 
MMT#07), spectra were obtained with two different 
MMT/Hectospec fiber configurations. In both cases, 
[O ill] A4363 was detected in each spectrum, confirming 
that these detections are robust. For these two galaxies, 
we combine the spectra for a higher S/N spectrum. In 
addition, more recent MMT/Hectospec observations for 
MMT#01, #03-#05, #07, #10, and #13 were taken in 
less ideal observing conditions. These observations also 
detected [O ill] A4363 in MMT#04, #07, and possibly 
#13. Since these spectra are less sensitive than those 
obtained in 2008, we do not combine them for the pur- 
pose of our analyses. 


3.1. Contamination from LINERs and AGN 

Because [O ill] A4363 is more likely to be detected in 
higher temperature gas, a common concern is whether 
these galaxies harbor low-ionization nuclear emitting re- 
gions (LINERs; Heckman 1980), where the gas may be 
shock-heated. To determine if any of our galaxies are 
LINERs, the preferred method is to use the [O I] A6300 
emission line. Unfortunately, this is redshifted out of 
our spectral coverage for the majority of our galaxies. 
Instead, we use a variety of emission-line flux ratios, in- 
cluding [O 11 ] /[O hi], [O iii]/H/3, and [O n]/[Ne in] A3869, 
to determine if any of our galaxies could be a LINER. 
Four galaxies (MMT#02, #03, #06, and #12) have 
[O 11 ] /[O III] ratios that are similar to or above unity 
(0.82-2.2), which is a cautionary LINER flag. Upon com- 
paring our emission-line fluxes to SDSS DR7 LINERs, we 
find that only MMT#03 is arguably a LINER. The three 
remaining galaxies have too low (high) of an [O ll ]/[0 ill] 
([O m]/H/3) ratio by at least 0.4 dex. Further indepen- 
dent evidence supporting the idea that MMT#03 is a 
LINER is its stellar mass. As we will later show (Sec- 
tion 4.4), it is our most massive galaxy. In general, LIN- 
ERs are primarily found in more massive galaxies (e.g., 

13 There are two galaxies (Keck#2, #4) that we include in our 
sample for various reasons discussed in Section 4.1 and Table 3. 


94% of SDSS DR7 LINERs are above a stellar mass of 
10 10 Mq). 

Another possibility to consider is that some of 
our galaxies might harbor an AGN. Supporting evi- 
dence would include very high ionization lines (e.g., 
[Ne V] A3425, He II A4686), although these lines have also 
been seen in some local blue compact dwarf galaxies (Izo- 
tov et al. 2012b). A search for [Ne V] and He II yielded 
non-detections, arguing that our sample is free of Seyfert 
galaxies. In addition, our two galaxies (MMT#10 and 
#11) with coverage of [N 11 ] A6583 have [N n]/H a ra- 
tios of 0.02 and 0.04, which places them in the star- 
forming region of the “BPT” diagram (Baldwin et al. 
1981). While [N 11 ] is not available for the remaining 18 
galaxies, we are able to use an analogous diagnostic tool 
called the “Mass- Extinction” diagram (“MEx”; Juneau 
et al. 2011). We illustrate in Figure 4 our sample against 
star-forming galaxies selected using the BPT diagram 14 
from the SDSS MPA-JHU DR7 sample 15 . We find that 
all of our galaxies lie in the star-forming domain. Com- 
pared to local galaxies of similar stellar masses, our ELGs 
all show a higher [O III] /H/3 flux ratio by « 0.2-0. 5 dex. 
We find that for more typical ELGs (i.e., those that lack 
[O in] A4363 detections) in our spectroscopic sample, the 
[O ill] /H/3 ratios are similarly higher than in SDSS. This 
holds for more massive galaxies with M* ~ 10 10 Mq. We 
overlay the range of the full sample as a purple shaded 
region in Figure 4. We later discuss this result for higher 
ionization in Section 5.4, and discuss selection effects as- 
sociated with our sample in Section 5.2. 

4. RESULTS 

In this section, we utilize our spectroscopic and pho- 
tometric data to estimate dust attenuation (§4.1), elec- 
tron temperature and gas-phase metallicity (§4.2), de- 
reddened SFRs (§4.3), stellar properties from SED mod- 
eling (§4.4), the local environment (§4.5), and the SFR 
surface density (§4.6). 


4.1. Dust Attenuation Correction from Balmer 
Decrements 


To correct the emission-line fluxes for dust attenua- 
tion, we use Balmer decrement measurements obtained 
from a combination of our spectroscopy and NB imag- 
ing. Since our ELGs possess high EWs, 12, 18, and 20 
galaxies have HJ, H 7 , and H/3 detected at >5<r, respec- 
tively. In addition, Ha measurements are available for 5 
galaxies. Two of them (MMT# 10 and MMT#11) are at 
lower redshifts, allowing us to use their spectroscopic Ha 
measurements. For the other 3 galaxies (MMT#04, #06, 
#11), the NB Ha flux is determined using the following 
equation: 


P A ATP /a,NB /a.BB J. 

Ha “ NB 1- (ANB/ABB)^' 


(6) 


where /a,nb and /a,bb are the flux density in erg s 1 
cm -2 A ' 1 for the narrow-band and broad-band, A’s are 
the respective FWHM of the filters (A z' = 956A), and 


14 We require at least 3 a detections of Ho;, [N il] , [O ill], and 
H/3, which yielded 274,613 galaxies. The sample was further lim- 
ited to 203,630 galaxies by excluding AGN and LINERs with the 
Kauffmann et al. (2003) selection. 

15 http : / / www . mpa-gar ching . mpg . de/SDSS/DR7 / . 
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TABLE 2 

Summary of Spectroscopic Sample 


ID 

Name 

Line Sel. 

R.A. 

Dec. 

^spec 

Obs. Dates (UT) 

^int 




[hr] 

[deg] 




[min.] 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


(7) 

(8) 

MMT#01 

NB816-140623 

[O in] 

13:25:16.87 

27:39:06.92 

0.6380 

2008 

Apr 14 

120 

MMT#02 

NB711-064628 

[O in] 

13:23:39.13 

27:32:52.71 

0.4327 

2008 

Mar 13 

120 

MMT#03 a 

NB973-104154 

Ha 

13:23:39.17 

27:31:47.34 

0.4809 

2008 

Mar 13 

120 

MMT#04 

NB704-088982.NB921-126525.IA679-112491 

[O III], Ha 

13:24:46.63 

27:34:56.98 

0.3933 

2008 

Mar 13, Apr 11 

240 

MMT#05 

IA679-031637 

[O m] 

13:25:03.37 

27:17:23.77 

0.3846 

2008 

Mar 13 

120 

MMT#06 

NB704-036405.NB921-063205 

[O III], Ha 

13:23:54.77 

27:20:12.53 

0.3995 

2008 

Mar 13 

120 

MMT#07 

NB704-049936-NB921-079428-IA679-062450 

[O III], Ha 

13:24:06.94 

27:24:01.76 

0.3896 

2008 

Mar 13, Apr 11 

240 

MMT#08 

NB816-081644 

[O III] 

13:23:42.97 

27:26:35.59 

0.6335 

2008 

Apr 10 

80 

MMT#09 

NB973-094500 

Ha 

13:23:49.80 

27:28:35.32 

0.4788 

2008 

Apr 10 

80 

MMT#10 

NB704-009999 

Ha 

13:23:56.39 

27:13:32.98 

0.0683 

2008 

Apr 10 

80 

MMT#11 

IA598-079010 

[O in] 

13:24:13.64 

27:25:09.27 

0.1752 

2008 

Apr 10 

80 

MMT#12 

NB816-112403 

[O III] 

13:25:21.78 

27:33:15.69 

0.6405 

2008 

Apr 10 

80 

MMT#13 

NB711 102472.NB973-156739 

H/3, Ha 

13:24:28.88 

27:45:51.88 

0.4696 

2008 

Apr 11 

120 

MMT#14 

NB711-077774_NB973-125003 

H/3, Ha 

13:25:22.94 

27:37:40.33 

0.4644 

2008 

Apr 11 

120 

Keck# 1 

NB711-049857.IA679-079866.NB921-995851 

[O in] 

13:25:11.94 

27:27:31.20 

0.8390 

2004 

Apr 23 

120 

Keck# 2 

NB816-070113 

[O in] 

13:24:34.91 

27:24:10.20 

0.6230 

2004 

Apr 23 

118 

Keck# 3 

IA679-077341 

[O 11 ] 

13:23:53.54 

27:27:13.01 

0.7906 

2004 

Apr 23 

118 

Keck# 4 

NB 704-087569 

[O II] 

13:23:43.51 

27:34:20.98 

0.8829 

2008 

May 1 

130 

Keck# 5 

NB921-078003 

[O III] 

13:24:43.66 

27:23:34.86 

0.8353 

2009 

Apr 27 

180 

Keck# 6 

NB704-060432 

[O m] b 

13:24:58.62 

27:26:40.47 

0.8237 

2009 

Apr 27 

180 


a This source is likely a LINER (see Section 3.1). 

^This galaxy was previously targeted for its NB921 excess flux due to [O III]. It was not identified as a NB921 emitter in the re-selection because of 
its faintness. The NB704 excess is due to [Ne III] A3869. 


TABLE 3 

Emission-line Measurements and Flux Ratios 


ID 

[O ii] 

H/3 

[O hi] A4959 

[O hi] A5007 

[O hi] A4363 

S/N(4363) 

7?23 

032 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

MMT#01 

11.05i0.59 

9.55=1=0.78 

19.12il.04 

54.30i0.70 

2.20i0.50 

4.44 

8-851“* 

6 64 +0 ‘ 44 

MMT#02 

18.55i0.35 

5.93=1=0.27 

6.73i0.27 

19.41i0.28 

0.64i0.21 

3.06 

7-54l“i* 

1 41+0.05 
i41 -0.05 

MMT#03 a 

20.94±0.46 

5.79=1=0.27 

3.13i0.39 

9.51i0.31 

0.84i0.21 

3.95 

5.80tg;« 

0 60 +0 ‘ 03 
U.DU _ 0 03 

MMT#04 

50.23±0.28 

22.17=1=0.17 

32.13i0.26 

99.65i0.24 

1.48i0.22 

6.86 

o o-| +0.05 
°- z± -0.06 

9 ^9+0.02 

Z.UZ-O 02 

MMT#05 

3.84i0.36 

7.09=1=0.25 

15.15i0.24 

43.97i0.26 

1.30i0.25 

5.22 

o qq+0.40 
o.oO_q 32 

15.3811:8® 

MMT#06 

21.76±0.34 

6.31=1=0.34 

5.66i0.29 

16.34i0.26 

0.84i0.27 

3.10 

6-941°;" 

l-Ollara 

MMT#07 

20.15±0.28 

21. 06=1=0. 20 

35.97i0.24 

120.77i0.24 

2.28i0.26 

8.91 

8-40l°;i° 

7 7Q+0-11 
' • ' 0.08 

MMT#08 

28.23i0.71 

11. 21=1=1. 27 

10.47il.04 

32.57i0.89 

1.85i0.56 

3.32 

6-36l°;** 

1 c 9 +0.08 

MMT#09 

16.88±0.82 

6.81=1=0.62 

10.37il.12 

28.96il.17 

1.45i0.41 

3.53 

8 - 261 °+ 

9 oq+0.05 
z.oo_ 0 10 

MMT#10 

22.70il.06 

18.55i0.42 

30.93i0.45 

92.80i0.52 

2.22i0.62 

3.58 

7-891°+ 

"5 4^+0. 21 
0-40_o.25 

MMT#11 

43.34i0.75 

22.30i0.57 

35.26i0.60 

112.99i0.63 

2.08i0.45 

4.67 

8-5918:12 

3 49+ 0 06 
0.07 

MMT#12 

24.18i0.68 

6.30i0.93 

11.29i0.77 

29.59i0.73 

1.46i0.48 

3.04 

10.32li;|| 

1 69 - *" 0-07 

l.DU_o 07 

MMT#13 

6.33±0.54 

6.55i0.26 

10.94i0.39 

28.50i0.45 

1.07i0.31 

3.46 

6 99 + 0 33 

u.»»_o 34 

c 9q+0.66 

0.63 

MMT#14 

102.13i0.54 

57.20i0.34 

87.06i0.61 

269.12i0.67 

2.14i0.28 

7.55 

8-011^6 

3 49 + 0 - 02 
0.02 

Keck#l 

4.61i0.14 

7.72i0.09 

15.35i0.09 

45.52i0.08 

0.74i0.07 

11.13 

8-49l°;J° 

13 21 + 0 - 46 

lo.Zl_o 45 

Keck# 2 

3.55±0.15 

2.54i0.05 b 

4.31i0.05 

13.40i0.05 

0.64i0.10 

6.22 

8-3618:11 

3 no +0 - 30 

O.UU_o 29 

Keck#3 

27.89i0.16 

12.60i0.ll 

14.78i0.07 

44.94i0.08 

0.30i0.06 

4.64 

6-95l°;°® 

2 14+0- 02 
z,± ^- 0.02 

Keck#4 

13.86i0.14 

7.37i0.11 

8.66i0.14 

26.85i0.14 c 

0.34i0.10 

3.44 

6.70l°;°3 

2 56+ 0 - 04 

z.ou—o.04 

Keck# 5 

1.86±0.07 

1.07i0.04 

1.63i0.04 

6.43i0.05 

0.55i0.16 

3.46 

9 30” 1 ” 0-48 

».ou_0 40 

4 04+O.I8 
0.18 

Keck #6 

<0.07 

0.92i0.03 

1.44i0.03 

4.00i0.03 

0.21i0.05 

4.27 

5 9fi “*" 0-20 
o.yo_ 0 16 

>80.70 


Note. — All emission-line fluxes are reported in units of 10 17 erg s 1 cm 2 with 68% confidence uncertainties. No dust attenuation 
corrections have been applied to these fluxes or flux ratios. 

“This source is likely a LINER (see Section 3.1). 

^The H/3 line fell inside a CCD gap. The reported flux here is derived from H 7 (see Section 4.1). 

c [0 III] A5007 was outside the spectral coverage; however, [O III] A4959 was detected at ~60cr. We therefore adopt a flux that is 3.1 
times the 4959A flux. 


Extremely Metal-Poor Emission-line Galaxies 


7 




0.6 

0.4 

0.2 

0.0 

- 0.2 


“I 1 | | I. I I I | I I I I 

Keck#2 6.2a 




I 1 I I I I, I I I I I I I I 








I I I I I I I I 1 


| i i i | | i, i i i | i i — r 

Keck#1 11.1a 






l l l l I l. I l l I l l l l 


i. i i i [ i r i i 

Keck#6 4.3a 




1 1 1 1 I ll 1 1 1 I 1 1 1 1 


4300 4350 4400 44504300 4350 4400 44504300 4350 4400 44504300 4350 4400 44504300 4350 4400 4450 


Rest Wavelength (Angstroms) 

Fig. 2. — Detections of [O ill] A4363 in 20 galaxies. The rest-frame spectra are shown in black, with vertical red dashed lines indicating 
the locations of H 7 A4340 and [O ill] A4363. OH sky-lines are indicated by blue vertical bands, where the strength of the sky-lines is denoted 
by their shade of blue (darker is stronger). The grey-shaded regions indicate the locally measured lcr noise, while the horizontal brown 
lines correspond to the median of the continua. 
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Fig. 3. — Emission-line luminosities and rest-frame EWs for our [O Ill]-A sample. All fluxes are observed, before correction for dust 
attenuation. Gray squares illustrate the SDSS DR7 emission-line sample with at least 5rr detections in all four lines. The density of sources 
is shown with a linear scale from white (low) to black (highest). Our only LINER (MMT#03) is denoted by the orange crosses. For 
comparative purposes, we show the EW and luminosity limits of our NB921 imaging (solid red lines) for z = 0.4 Ha (top left) and z = 0.84 
[O III] (bottom left) and NB704 imaging (solid blue lines) for z = 0.4 [O III], The MMT sample is labelled in blue while the Keck sample 
is indicated in green. 


/corr is the correction for the non-tophat shape of the 
NB filter when the redshifted Hct emission is in the fil- 
ter’s wing. Using filter throughputs and spectroscopic 
redshifts, we are able to determine f COII . We note that 
for four galaxies (MMT#03, #07, #13, and #14), the 
NB fluxes are not reliable for Balmer decrement determi- 
nations. In particular, MMT#03 suffers from significant 
[N ii] contamination in the NB921 filter as a LINER can- 
didate. For the other three sources, the Ha line falls in 
the wing of the NB921 or NB973 filter where precise filter 
response corrections cannot be made. 

A significant problem encountered with using Balmer 
decrements to determine dust attenuation is the under- 
lying stellar absorption. In three galaxies (MMT#03, 
MMT#04, and Keck#3), the spectra are of high S/N 
to sufficiently detect the continuum and fit and remove 
the stellar absorption or use IRAf’s splot command to 


re-measure the continuum from the absorption trough. 
For MMT#03 (Keck#3), we determine corrections of 
EW abs (H6) = 2.7A (1.6A), EW abs (H 7 ) = 1.9# (2.5#), 
and EW abs (H/3) = 1.6# (O.OA). While for MMT#04, we 
measure EW abs (Hi5) = 3.9A and EW abs (H 7 ) = 2.5A. 

For the remaining galaxies, the S/N of our spectra 
are not sufficient to model the stellar absorption. To 
correct these galaxies, we adopt EW abs (H(5) = 2A and 
EW abs (H 7 ) = lA. We assume no stellar absorption for 
H f3 and Ha, which is reasonable since the measured rest- 
frame emission-line EWs are significantly large (H/3: me- 
dian of 60A, average of 95A). 

With these corrections for stellar absorption, we illus- 
trate the Balmer decrements in Figure 5. In addition, we 
report in Table 4 the absorption-corrected Balmer decre- 
ments, and the determined or assumed EW for stellar 
absorption. 
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Fig. 4. — The MEx diagram for our [O ill]- A sample (circles) plot- 
ted against the SDSS DR7 emission-line sample (grey boxes). Each 
of our galaxies is labelled by a two-digit (MMT; blue) or one-digit 
(Keck; green) number. The stellar masses are determined from 
SED fitting, discussed in Section 4.4. The solid black lines defined 
by Juneau et al. (2011) separate star- forming (left), AGN (right), 
and composite (middle) galaxies. All of our galaxies are classified 
as star-forming based on the MEx classification. Our ELGs have 
[O III] /H/3 fluxes that are 0.2-0. 5 dex higher than the average for 
SDSS galaxies with the same stellar mass. We also show the SDF 
spectroscopic sample of ~200 z = 0.4-1 galaxies in purple with the 
shaded region representing ±lcr. This comparison indicates that 
ELGs at these redshifts have higher ionization in general, across 
2.5 dex in stellar mass. Our only LINER (MMT#03) is denoted 
by the orange cross. 

Under the assumption that the hydrogen nebular 
emission originates from an optically thick ionization- 
bounded H II region obeying Case B recombination, 
the intrinsic Balmer flux ratios are: (Hct/H/3)o = 2.86, 
(H 7 /H/3) 0 = 0.468, and (H<5/H/3) 0 = 0.259 for T e = 10 4 
K and electron density of n e = 100 cm' 3 (see Section 4.2 
for further discussion). We note that these values differ 
by less than 5% for T e = 2 x 10 4 K. Dust absorption 
alters these observed ratios: 


(Hn/H/3) obs = 

i Q-0.4E(B-y)[fe(Hn)-fe(H/3)] 

(Hn/H/?)„ 


( 7 ) 


where E(B-V) is the nebular color excess, and k( A) = 
A(X)/ E(B-V) is the reddening curve at A. The latter 
is dependent on the dust reddening “law.” We illustrate 
in Figure 5 the observed Balmer decrements under the 
Calzetti et al. (2000) (hereafter COO) dust reddening for- 
malism with fc(Ha) = 3.33, £;(H/3) = 4.60, fc(H 7 ) = 5.12, 
and fc(Hd) = 5.39. We find that our Balmer decrements 
are consistent with the COO dust reddening formula. For 
the remainder of our paper, all dust-corrected measure- 
ments adopt COO reddening. 

Our color excesses, which are tabulated in Table 4, 
are determined using either Ha/H j3 (#06, #09-#ll) or 
H 7 /H/3 (MMT#02-#05, #08, #12-#14, and Keck#l, 
#3, and #4). We do not use H<5/H/3 since it is not 
well measured for the majority of our galaxies. How- 
ever, the E(B-V) estimates from Hd/H/3 are consistent 
with H 7 /H /3, as indicated in Figure 5. In three galaxies 
(MMT#04, MMT#05, and Keck#4), the H 7 /H/3 ratios 
suggests negative reddening, so we adopt E(B-V) = 0.0. 


For four of our galaxies (MMT#01, #07, #5— #6), the 
dust reddening could not be determined from Balmer 
decrements. For these galaxies, we assume E(B~V) = 
0.24 ± 0.21 mag (A(Ha) ss 0.8 ± 0.7 mag), which is the 
average of our sample 16 , and is similar to typical redden- 
ing found for local galaxies (A(Ha) = 0. 8-1.1; Kennicutt 
1998a). This also agrees with what we found for ELGs 
at z « 0.4 (Ly et al. 2012a). 

For Keck#2, the H/3 line unfortunately fell at the edge 
of a CCD gap, so the H/3 flux is not fully measured. How- 
ever, the H 7 and H(5 lines are robustly detected (19 and 
9cr). Assuming Case B recombination and no reddening, 
the H 7 /H/3 and H(5/H/3 Balmer decrements yielded H/3 
fluxes of 2.54 x 10” 17 and 2.87 x 10” 17 erg s” 1 cm” 2 , 
respectively. This excellent agreement in predicted H/3 
fluxes suggests that very little reddening is present in this 
galaxy. We assume E(B-V) = 0.0. 


4.2. T e and Metallicity Determinations 

To determine the gas-phase metallicity for our galax- 
ies, we use the empirical relations of Izotov et al. (2006a). 
This follows the approach of most direct metallicity stud- 
ies. The first equation estimates T e ([ O III]) using the 
nebular-to-auroral [O ill] flux ratio: 


/ [O in] AA4959, 5007 3 
° S V [O hi] A4363 , 

where *3 = T e ([0 iii])/10 4 K, 
C T = (8.44 - 1.09*3 + 0.5 1\ 


1.432 

*3 


+ log Ct, 


0.08*1) 


1 + 0.0004a; 

1 + 0.044a; ’ 


(8) 

( 9 ) 


and x = 10” 4 n e *3 °' 5 . For the majority of our galax- 
ies, the [S Ii] AA6717,6732 doublet, an estimate for n e , 
is redshifted out of our spectral coverage. For one of 
our galaxies, MMT#10, both [S n] lines are weakly de- 
tected. The A6717/A6732 flux ratio of wl.l corresponds 
to n e ~ 400 cm” 3 (assuming T e = 10 4 K). In addi- 
tion, DEIMOS has the spectral resolution to separate 
[O n] AA3726,3729. The A3729/A3726 flux ratios vary 
between 0.87 and 1.35, which correspond to n e = 70-600 
cm” 3 . To determine T e , we assume n e = 100 cm” 3 . We 
note that Ct is only affected by n e in the high density 
regime (n e > 10 4 cm” 3 ), and therefore assuming n e = 
10, 100, or 10 3 cm” 3 yields nearly identical T e . 

We correct the nebular-to-auroral [O ill] flux ratio for 
dust attenuation using our dust attenuation prescriptions 
(see Section 4.1). Correcting for dust attenuation in- 
creases estimates of T e , since dust extincts shorter wave- 
lengths (e.g., 4363A) more than longer wavelengths (e.g., 
5007A). 

In addition, it has been known that T e determinations 
using (1) various direct metallicity prescriptions, includ- 
ing Izotov et al. (2006a); and (2) strong- line metallic- 
ity diagnostics do not yield consistent results (Kewley 
& Ellison 2008). These discrepancies have been recently 
explored, and it is believed that outdated effective col- 
lision strengths for the various 0 ++ excitation states, 
and a non-equilibrium distribution of the electron ener- 
gies result in overestimation of the electron temperature 
(Nicholls et al. 2012, 2013). Those authors provided ap- 
plicable corrections for each method. We adopt their 


16 The median E(B-V) is 0.25 mag. 
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Fig. 5. — Balmer decrements for our [O ill] A4363 galaxies. The MMT and Keck samples are shown as blue circles and green squares, 
respectively. Red crosses indicate that the H<5 fluxes are unreliable since it fell on an OH sky-line for MMT#07 and Keck#4. Vertical 
blue lines indicate non-detections in H5 for MMT#02 and MMT#09. Finally, red circles and the solid red curves show the effects on the 
Balmer decrements with increasing dust reddening following COO. The values reported in red are A(Ha). 


TABLE 4 

Balmer Decrement Measurements and Derived Reddening 


ID 

Ha/H /3 

H7/H/3 

EW(H< 5 ) 

EW(H 7 ) 

EW(H/ 3 ) 

Source 

E ( B - V ) 

(1) 

(2) 

( 3 ) 

( 4 ) 

( 5 ) 

(6) 

( 7 ) 

(8) 

MMT #01 


a 

2.0 

1.0 

0.0 


0 24“*" 0,21 
u ‘ z ^-0.21 

MMT #02 


0 41 + 0-03 

u - 4 i - 0.04 

2.0 

1.0 

0.0 

H 7 /H£ 

0 o^”*" 0 ’ 20 
u.zo_022 

MMT #03 

b 

q og+ 0.04 
u.o»_o 04 

2.7 

1.9 

1.6 

h 7 /h^ 

U.OO-o 24 

MMT #04 

3.98±8:SS- 

0 47 + 0 - 01 

u - 4 '-o.oi 

3.9 

2.5 

0.0 

h 7 /h^ 

o.ool°;°® d 

MMT #05 


0 49 -1 " 0,05 
u, ^ y - 0.04 

2.0 

1.0 

0.0 

H 7 /H /3 

0.00l° : ^ d 

MMT #06 

4-5i±g:Ii 

O q' 7 + 0-06 

u, °' - 0.05 

2.0 

1.0 

0.0 

Ha/H /3 

0 39 +0 04 

u.o »_0 03 

MMT #07 

C 

0 . 48 l °;°2 

2.0 

1.0 

0.0 


0 24”*” 0,21 
U.Z4-0 21 

MMT #08 


0 - 38 l °;°5 

2.0 

1.0 

0.0 

h 7 /h^ 

0 . 42 l°;f 5 

MMT #09 

4 45 + 0 - 45 
^°- 0.44 

0 36 +010 
u.oo_ 0 1Q 

2.0 

1.0 

4.4 

Ha/H p 

o. 38 i°;°° 

MMT #10 

3 . 26 lg?f 

0 47 + 0 - 04 
- 0.04 

2.0 

1.0 

0.0 

Ha/H /3 

n 1 1+0.02 
U,11 -0.02 

MMT #11 

3 59 ”*”°' 12 
o.o»_ 0 12 

0 41+ 0 03 
0.02 

2.0 

1.0 

0.0 

Ha/H /3 

0 19 - *" 003 
0.03 

MMT #12 


0 47 + u - u7 

u - J '-o.n 

2.0 

1.0 

0.0 

H 7 /H^ 

U - O 1 - 0.46 

MMT #13 

c 

0 43 + U U5 
d. z ±o_ 0 05 

2.0 

1.0 

0.0 

h 7 /h^ 

0 17+0.20 
u, 1 / - 0.24 

MMT #14 

d 

0.46l°;°oi 

2.0 

1.0 

0.0 

h 7 /h^ 

0 03+0 03 
U.OO-o 03 

Keck#l 


0 . 38 l°;°i 

2.0 

1.0 

0.0 

h 7 /h /3 

0 41 + 0 - 07 
u - 41 - 0.07 

Keck# 2 


a 

2.0 

1.0 

0.0 


0.00 e 

Keck #3 


0 315+ 001 
D.OO_ 0> 01 

1.6 

2.6 

0.0 

H 7 /H /3 

0.621°;°* 

Keck #4 


0 53+ 0 01 

U.OO -001 

2.0 

1.0 

0.0 

h 7 /h^ 

0.00l°;°^ d 

Keck #5 


a 

2.0 

1.0 

0.0 


0 24”*” 0,21 
u ‘ z ^-0.21 

Keck#6 


_ f 

2.0 

1.0 

0.0 


0 24“*" 0,21 
u ‘ z ^-0.21 


Note. — Uncertainties are reported at the 68% CL. Stellar absorption line corrections are reported 
in Cols. (4)-(6) 

a H/3 is affected by a weak OH sky-line. 

^The NB921 flux is contaminated by significant [N II] emission because this source is a LINER. 
c The Ha emission falls at the very edge of NB921 or NB973 filter, and therefore extinction deter- 
minations are unreliable using Ho:/H/3. 

^Negative reddening was seen. We adopt E(B—V ) = 0.0. 
e The H7 and H<5 fluxes suggests no reddening 
/h 7 is affected by an OH sky-line. 
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corrections, which correspond roughly to a 5% reduction 
for T e estimates using the Izotov et al. (2006a) approach. 

For the very strong ELGs we are studying, our [O ill] 
measurements have a very large dynamic range. The 
strongest [O ill] A4363 line is as much as 0.089 of the 
[O in] flux (MMT#03), while the weakest is 0.007 
(Keck#3). We find that the average and median 
A4363/A5007 flux ratio for our sample are both 0.037. 
The derived electron temperatures for our galaxy sam- 
ple, also have a wide range, from 10 4 K to 4.5 x 10 4 
K. Among our 19 galaxies, 11 have T e estimates of 10 4 - 
2.1 x 10 4 K, which are similar to those measured in local 
galaxies 1 '. The large tail toward high T e yields an av- 
erage log(T e /K) that is higher than the average of local 
galaxies by 0.15 dex. 

Other higher redshift studies have measured high T e 
in star-forming galaxies. For example, Yuan & Kewley 
(2009) detected [O ill] A4363 in a lensed z = 1.7 galaxy 
and measured T e ss 2.35 x 10 4 K, and Hu et al. (2009) 
have also measured large [O ill] A4363/A5007 flux ratios 
in a few ELGs. The higher values of T e are not well- 
determined, and have a large tail (0.2-0. 3 dex) toward 
lower temperatures. This is likely the case for other pub- 
lished measurements, as Yuan & Kewley (2009) yielded 
a 3cr detection of [O ill] A4363 and the measurement un- 
certainties of Hu et al. (2009) result in a significant tail 
toward low T e . 

MMT#03 has the most extreme T e value, but as we 
mentioned in Section 3.1, it is a LINER with shock- 
heated, rather than photo-ionized gas producing most 
of its emission-line spectrum. Other galaxies with high 
T e (> 2.5 x 10 4 K) are MMT#06, #08, #09, #12, 
Keck#5, and #6. To examine if these temperatures 
are reasonable, we generate low-metallicity models with 
CLOUDY (Ferland et al. 1998). In our modeling, we in- 
put Starburst99 (Leitherer et al. 1999) SEDs that adopt 
Geneva stellar evolution models with five stellar metallic- 
ities between Z I Zq = 0.05 and 2.0. Coupling the stellar 
and gas metallicity, we find that the highest plausible 
T e ([0 III]) temperature is «25,000 K for Z / Zq = 0.02- 
0.03. 18 Given this upper limit, we fix the dust-corrected 
T e ([0 ill]) to not exceed this temperature for these six 
galaxies. We note that without this T e offset, the result- 
ing oxygen metallicity will be lower by 0. 1-0.4 dex for 
these six galaxies. 

The T e measurements for all 20 objects are shown in 
Figure 6. Throughout this paper, we generate realiza- 
tions of the emission-line fluxes to construct probability 
distribution functions (PDFs) for all observed and de- 
rived measurements. This is critical, as the distributions 
are non-Gaussian in the domain of low S/N (<5), and 
it allows us to propagate our measurement uncertainties, 
including E(B-V). 

To determine the ionic abundances of oxygen, we use 
two emission-line flux ratios, [O n] AA3726,3729/H/3 and 

17 90% of Nagao et al. (2006) ’s sample has T e = 9,400-20,000 
K. 

18 An extrapolation is made from Z/Zq = 0.05. 


[O III] AA4959,5007/H/3: 

— 0.41ogt 2 ~ 0.034t 2 + log(l + 1.35a:) 
-0.551ogf 3 -0.014*3. 


Here, t 2 refers to the singly-ionized oxygen electron tem- 
perature, T e ([ O n] ) . For our metallicity estimation, we 
adopt a standard two-zone temperature model with t 2 = 
T e ([0 ii])/10 4 K = —0.577 + t 3 (2.065 - 0.498t 3 ) (Izotov 
et al. 2006a). In computing 0 + /H + , we also correct the 
[O n]/H/3 ratio for dust attenuation. We do not cor- 
rect 0 ++ /H + since the effects are negligible (e.g., the 
[O iii]/H/3 ratio changes by less than 0.02 dex for A{ Ha) 
= 1 mag). 

Since the most abundant ions of oxygen in H II re- 
gions are 0+ and 0 ++ , we can determine the oxygen 
abundances as: 


O 

H 



( 12 ) 


In Table 5, we provide observed and de-reddened 
flux ratios, estimates of T e ([ O ill]), log(0 + /H + ), 
log(0 ++ /H + ), and 12 + log(0/H) for our sample, and 
95% confidence uncertainties (i.e. , two standard de- 
viations). In addition, we illustrate our PDFs of 
12 + log(0/H) in Figure 7. Our three most metal- 
poor galaxies are MMT#13, Keck#2 and #6 with 
12 + log(0/H) = 7.53±8;11, 7.54+o.22> and 7.24t°; 49 , re- 
spectively. For our sample of 19 star-forming galaxies, 4 
of them can be classified as an XMPG as their metallicity 
is below 12 + log(0/H) = 7.65. 


4.3. Dust- Corrected Star Formation Rates 

In addition to gas-phase metallicity determinations, 
our extensive data allow us to determine dust-corrected 
SFRs using the Hydrogen recombination lines, which are 
sensitive to the shortest timescale of star formation, <10 
Myr (i.e., an instantaneous SFR estimate). 

Assuming a Chabrier (2003) (hereafter Chabrier) ini- 
tial mass function (IMF) with minimum and maximum 
masses of 0.1 and 100 Mq, and solar metallicity (Ken- 
nicutt 1998a) 19 , the SFR can be determined from the 
dust-corrected Ha luminosity as: 

SFR(M 0 yr _1 ) = 4.4 x 10" 42 L(erg s -1 ). (13) 

Similarly, SFRs can be determined from H/3 by first cor- 
recting for dust reddening (A(H/3) = 4.60 E(B~V)) 1 and 
then adopting the intrinsic Case B flux ratio, (Ha/H/3)o 
= 2 . 86 . 

Our SFR estimates are summarized in Table 6 and 
illustrated in Figure 8. We find that our galaxies have 
dust-corrected SFRs of 0.04-64 Mq yr _1 with an average 
(median) of 8.1 (2.3) Mq yr _1 . 

4.4. Spectral Energy Distributions and Estimated Stellar 
Population Properties 


19 We assume a factor of 1.8 between the integrated masses for 
the Chabrier and Salpeter (1955) (hereafter Salpeter) IMFs. 
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Fig. 6. — Probability distribution functions for the electron temperature determined from the [O III] A4363 line. The peak values and 
±95% confidences are reported in each panel. T e ’s are derived from the observed flux ratios (hashed histograms) and the dust-corrected 
flux ratios (grey histograms). 


One significant advantage of studying low-mass galax- 
ies in the SDF is the ultra-deep imaging in twenty- 
four bands from the UV to the IR, which allows us to 
characterize their stellar properties. The SDF is im- 
aged with: (1) GALEX (Martin et al. 2005) in both the 
FUV and NUV bands 20 ; (2) KPNO’s Mayall using MO- 
SAIC in U ; (3) Subaru with Suprime-Cam in 14 bands 
(BVRci 1 z' ZbZ r , and the five NB and two IA filters as 
mentioned previously); (4) KPNO’s Mayall using NEW- 
FIRM (Probst et al. 2008) in J and H, (5) UKIRT us- 
ing WFCAM in K, and (6) Spitzer in the four IRAC 
bands (3.6/rm, 4.5/rm, 5.8^tm, and 8.0/rm). As discussed 
in Ly et al. (2011a), a more complete photometric cata- 
log from SExtractor was constructed from an ultra-deep 
mosaic that consisted of stacked optical and near-IR im- 
ages with a normalization that takes into account the 
rms (i.e. , y 2 stacking). 

Upon cross-referencing our 20 galaxies with the com- 
plete photometric sample, all but one of them have di- 
rect matches. The missing source (MMT#05) was con- 

20 Details on the GALEX imaging are available in Ly et al. 
(2009) and Ly et al. (2011a). 


fused with a nearby bright extended galaxy. We exclude 
MMT^05 in our stellar population analyses. For the re- 
maining 19 galaxies, aperture photometry measurements 
in 15 bands from the FUV to K are used to construct 
their complete broad-band SEDs. There are two modifi- 
cations that we make to these measurements. 

First, since all our galaxies are virtually point sources 
for GALEX , we obtain more accurate photometry in the 
FUV and NUV bands by PSF-fitting with daophot 
(Stetson 1987). Our examination of the residuals for 
each source suggests that the PSF-fitting is extremely 
successful, and that these galaxies are not affected by 
contamination of nearby sources. Among the Keck sam- 
ple, we note that only Keck#2 is detected in the FUV. 
The lack of FUV detections is un-surprising, since the 
Lyman continuum break occurs redward of the filter for 
the remaining five Keck sources. Also, Keck#6 is not de- 
tected in the NUV because of its intrinsic faintness. Be- 
cause the MMT [O ill] A4363 galaxies are brighter with 
z < 0.65, they are all detected in the NUV and FUV 
bands, allowing for robust UV SFR determinations. 

Second, since our galaxies all have very high emission- 
line EWs, we correct the broad-band photometry for the 
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Fig. 7. — Probability distribution functions for the gas-phase metallicity, 12 + log(0/H). The peak values and ±95% confidences are 
reported in each panel. We illustrate the observed and de-reddened determinations by the hashed and grey filled histograms, respectively. 
The upper limit of the [O II] flux is adopted for calculating the metallicity for Keck 41 6. 


contribution from nebular emission lines using emission- 
line measurements from our spectroscopy and NB imag- 
ing. Here we generate a spectrum for each galaxy with 
zero continuum and emission lines located at the red- 
shifted wavelengths for [O II], [Ne III], H/3, [O ill], Ha, 
and higher order Balmer lines. These spectra are then 
convolved with the filter bandpasses to determine excess 
fluxes, and then are removed from the broad-band pho- 
tometry. The correction for emission lines to SEDs has 
been recently explored in other higher redshift studies 
where high-sSFR galaxies are more prominent at 2 > 1 
(Atek et al. 2011; Pirzkal et al. 2013). For two galaxies 
(MMT#10 and MMT#11), the redshifts are low enough 
to have spectroscopic detections of Ha. For a subset of 
MMT galaxies (#03, #04, #06, #07, #09, #13, and 
#14), we use the NB921 or NB973 excess fluxes for Ha 
flux estimates. However, for galaxies at z > 0.5 and 
0.43 < z < 0.47, Ha is redshifted into the near-IR 
or not available from optical spectroscopy. To correct 
Ha in these galaxies (Keck#l-#6, MMT#1, #2, #8, 
and #12), we use spectroscopic information for H/3, but 
with the assumption that Ha is three times stronger. Of 
course, higher dust reddening will yield stronger Ha cor- 


rections, so we are adopting a minimum correction of 
A(Ha) = 0.14 mag. 

Both our observed SEDs and those corrected for nebu- 
lar emission lines are tabulated in Table 7, and represen- 
tative examples are shown in Figure 9. The black circles 
show the original broad-band photometry, while the blue 
circles show the corrected continuum fluxes after removal 
of emission lines. These three examples illustrate how 
the SEDs change when the emission-line corrections are 
small, average, and among the largest. 

Both the original and corrected SEDs were fitted with 
stellar population synthesis models (Bruzual & Char- 
iot 2003) with the Fitting and Assessment of Synthetic 
Templates (FAST; Kriek et al. 2009) code. We use 
exponentially-declining star formation histories (i.e., r 
models) similar to previous fitting by us and many other 
groups (Ly et al. 2011a). We have chosen this star- 
formation history (SFH) for its simplicity since broad- 
band data are generally unable to distinguish against 
more complicated SFHs (e.g., a constant SFR with a 
recent burst, which may be more representable for our 
galaxies). As we later find, these fits adopting an ex- 
ponentially declining SFH are consistent with the data 
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Fig. 8. — The SFR versus stellar mass for our [O Ill]-A sample. The stellar masses are obtained from SED fitting (Section 4.4) where we 
adopt a Chabrier IMF. The SFRs are determined from either Ha (filled squares) or H/3 (open squares), which are sensitive to a timescale 
of <10 Myr. The observed SFRs are shown as grey open squares with the degree of dust attenuation corrections indicated by the solid grey 
vertical lines. These dust corrections are discussed further in Section 4.1. The average (median) sSFR of our galaxies is 10~ 8 yr" 1 (8x 10 _J 
yr — 1 ). Each of our galaxies are labelled by a two-digit (MMT; blue) or one-digit (Keck; green) number. Our only LINER (MMT#03) is 
denoted by the orange cross. Overlaid in red crosses are the low-luminosity local sample of Berg et al. (2012), and the greyscale demonstrates 
the SDSS sample used by Mannucci et al. (2010). The “main sequence” relation of Noeske et al. (2007) at z = 0.45—0.7 and z = 0.7—0.85 
are illustrated in orange and brown, respectively. For direct comparisons, we have corrected Noeske et al. (2007) and Berg et al. (2012) 
measurements to a Chabrier IMF. 


with nearly unity y 2 values. Also, the primary purpose 
of our SED modeling is to obtain stellar masses of our 
galaxies. This is best traced from the rest-frame optical 
light from the older stellar population. Thus the inclu- 
sion of a recent burst does not significantly alter stellar 
mass estimates. 

In these models, we continue to adopt COO reddening, 
and inter-galactic medium attenuation following Madau 
(1995). The only differences are that we adopt a Chabrier 
IMF, allow for an extremely short burst model, r = 10 7 
yr, and use stellar atmospheres with one-fifth (Z=0.004) 
solar abundances. The latter is motivated by the fact 
that our galaxies are extremely metal-poor. In addition 
we considered solar abundance (Z= 0.02) for complete- 
ness. We find that for the majority of our galaxies, the 
^=0.004 models yield similar results to Z= 0.02 models. 
For consistency with our metallicity determinations, we 


adopt the SED fitting results with Z= 0.004 models. 

We find that with the emission-line corrections, the 
SED fits are improved with lower y 2 values by a factor 
1.1-10.4 in the Z= 0.004 models with an average (me- 
dian) improvement of 3.3 (2.1). In general, the results 
of fitting the SEDs with emission-line corrections yielded 
lower stellar masses by 0.1 dex and younger ages by 0.1 
dex on average. However, these differences vary from one 
source to another. In particular, our analyses also reveal 
a moderate correlation between [O III] EW and the dif- 
ference in stellar mass determinations, which others have 
reported (Atek et al. 2011). 

Our SED-fitting results are summarized in Table 6. We 
find that these galaxies typically are low-mass systems 
(median of 3 x 10 8 M@, average of 2 x 10 8 Mq). However, 
there is dispersion in the mass distribution (extending 
from 5x10® to 3x 10 9 Mq). The estimated light-weighted 
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Fig. 9. — The SED and the best-fitting stellar population results from modeling the SED in three galaxies: Keck#2, MMT#01, and #04. 
The observed fluxes are shown by the black circles while the blue circles illustrate the corrections for nebular emission lines. Errorbars 
along the ic-axis demonstrate the FWHM of the filters while those along the y-axis are 1 cr. The best x 2 fits to the original and emission-line 
corrected SEDs are shown in red and blue, respectively. In fitting the SED, we exclude the IRAC photometry, but include them in this 
figure for comparison. The emission-line corrected fitting results are reported on the right. Grey vertical lines indicate the locations of the 
Balmer and Lyman continuum breaks, while red, purple, and orange vertical lines show the location of various redshifted emission lines. 
When correcting for emission lines, the % 2 values are reduced by a factor of 3.5—10.3 in these three cases. In MMT#04, we find that the 
stellar mass and age are overestimated by 0.3 dex and 0.7 dex, respectively. 


stellar ages of log(t age /yr) = 7. 1-8. 7 dex (average: 7.75 
dex) suggest that these galaxies have undergone most of 
their star formation in the recent past. 

We compare the de-reddened SFRs estimated from 
SED fitting against the Balmer-derived de-reddened 
SFRs and find good agreement, albeit with large uncer- 
tainties. We also compare the dust reddening determina- 
tions from SED fitting and Balmer decrements and find a 
weak correlation. With large uncertainties in determin- 
ing reddening with both methods, it is not possible for us 
to investigate any difference between nebular and stellar 
reddening. The weak correlation is consistent with E{B- 
V) gas w 1.5E(B-V) a tai- This result is similar to what we 
found for a larger sample of ELGs in Ly et al. (2012a). 

Finally, we combine our dust-corrected SFRs (Sec- 
tion 4.3) and stellar mass determinations to illustrate 
our galaxies on the SFR-M* plane in Figure 8. Relative 
to their stellar mass, we find that our ELGs are all un- 
dergoing strong star formation. In particular, the sSFRs 
are between 10~ 9 yr -1 and 7 x 10~‘ yr -1 with an aver- 
age of 10” 8 yr -1 . This is 2. 0-3.0 dex higher than the so 
called “main sequence” of star formation for z ~ 0 SDSS 
galaxies (Mannucci et al. 2010) and low-luminosity local 
galaxies (Berg et al. 2012). Extrapolating the main se- 
quence of Noeske et al. (2007) toward lower stellar mass, 
we find that the sSFR of our ELGs are «0. 5-1.0 dex 


higher than “typical” galaxies at z ~ 0.45-0.85, but they 
are consistent with Noeske et al. (2007)’s sample at the 
la level. We emphasize that the inverse of the sSFRs, is 
consistent with the stellar ages derived from SED fitting. 

With the strong nebular emission lines seen in our 
galaxies, the nebular continuum from free-free, free- 
bound, and two-photon emissions can be important (see 
e.g., Izotov et al. 2011). To estimate the contribution to- 
ward the total light, at rest-frame optical wavelength, we 
generate spectral synthesis models from Starburst99. In 
our models we assume a constant SFH, a Kroupa (2001) 
IMF 21 , Geneva stellar evolutionary models, and Z / Zq = 
0.05. For EWo(H/3) = 45A (the median of our sample), 
we estimate that the nebular continuum is responsible 
for <10% of the total optical/IR light, and thus the stel- 
lar masses (sSFRs) that we have reported are <0.05 dex 
overestimated (underestimated). We note that this cor- 
rection can be as large as 20-50% (0. 1-0.3 dex) in galax- 
ies (Keck#l, MMT#07) with EW 0 (H/3) of «250-400A. 
However, these corrections for only two galaxies do not 
significantly alter our MZR. 

4.5. Nearby Emission-line Galaxy Companions 

21 This is similar to a Chabrier IMF. 
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Utilizing our deep NB/IA images, we constructed 
postage stamps that contain only the emission-line ex- 
cess flux by differencing the NB/IA images with adja- 
cent broad-band images (as previously discussed in Sec- 
tion 2.1). For proper removal of the continuum, we nor- 
malize the images by their respective zeropoints. We find 
that the majority (14 of 20) of our sample has between 1 
and 4 NB/IA excess emitters within a projected radius 
of sslOO kpc. The closest separation is 12 kpc with an 
average projected distance of 48 kpc. 

To ensure that these nearby NB /IA excess sources are 
real, we conducted two checks. First, we examine if 
any of these nearby galaxies are also classified as an 
NB/IA excess emitter in our source catalogs. We find 
that more than half of them (16 of 34) are NB/IA excess 
emitters. Upon further inspection, we find that those 
that are missed did not satisfy our NB/IA excess se- 
lections because they have low emission-line EW or are 
faint, but many could also lie at the same redshift as 
our [O ill] A4363 galaxies. Second, some of our galaxies 
are identified as an excess emitter in two or more filters 
(e.g., NB704 and NB921). This is due to multiple emis- 
sion lines falling in these filters. We find that in these 
four galaxies, six of the nearby excess emitters are also 
detected in two or three filters. This result strongly sug- 
gests that these nearby NB /IA excess emitters surround- 
ing our [O ill] A4363-selected galaxies are physically real, 
and at a similar redshift. 

While our sample size is small, the presence of nearby 
(<100 kpc) companions in 70% of our galaxies suggests 
that this is unlikely to be a coincidence. To confirm this, 
we examine how frequent nearby emission-line compan- 
ions are using our sample of 401 z = 0.4 Ha emitters, 
which has been previously discussed in Ly et al. (2012a). 
We determine the projected distance to the nearest Ha 
emitter, and find that the average of this distance is ~350 
kpc (median is «280 kpc). Furthermore we find only 
49 of 401 (12%) Ha emitters have a companion that is 
within 100 kpc. This low percentage suggests that Ha 
emitters, while clustered, are infrequently found with a 
nearby satellite. An identical analysis with 715 z = 0.84 
[O in] emitters yielded similar numbers. 

Three plausible explanations for the significant excess 
of nearby companions are that galaxy-galaxy interactions 
can: (1) stimulate a “starburst”, (2) tidally strip out 
metal-rich gas into the CGM, or (3) cause accretion of 
metal-poor gas from the outskirts of the galaxy. 

In the first case, the intense star formation can trig- 
ger strong winds that can entrain metal-rich gas into the 
CGM. If so, this would explain why these galaxies are 
identified to be extremely metal-deficient. Our measure- 
ments of SFR and masses (see Figure 8) suggest that the 
current SFR could form most of the stellar population in 
~100 Myr (average), supporting this hypothesis. Since 
our [O ill] A4363 sample is inhomogenous (i.e. , selection 
in different NB /IA filters) , a more detailed study is war- 
ranted from the complete NB/IA sample. Forthcoming 
work will examine further these nearby companions and 
any correlation seen with the properties of the central 
galaxies. 

In the case of “tidal disruption,” one would expect to 
see evidence of gas stripping, in the form of tidal tails 
and asymmetric morphologies. An examination of our 
NB/IA-excess images reveal no tidal tails or any “bridge” 


between our targeted sources and a nearby companion. 
We do notice that the light distribution is asymmetric 
in half of our galaxies; however, the seeing-limited res- 
olution of our data makes it difficult to interpret such 
results. Higher resolution imaging with Hubble Space 
Telescope (HST) is needed to quantitatively measure the 
morphological properties of these galaxies and detect any 
tidal tail features. 

One of our galaxies (Keck#05) fortuitously has 
HST/WFC3 broad-band near-IR imaging (Jiang et al. 
2013). Keck#05 is in fact two galaxies that are separated 
along the NE-SW direction by 0.9" (projected distance 
of 6.9 kpc at z = 0.84) with the nebular emission coming 
from the SW component. Our examination of the HST 
imaging reveal that the NE galaxy is diffuse along the 
P.A. toward the SW component. More careful analysis 
is needed to understand Keck#05 and its nearby com- 
panion, and as we emphasized, HST imaging for a larger 
fraction of our [O ill]- A sample is required to examine the 
possibility of tidal disruption for the case of low metallic- 
ities in our galaxies. We note however that this scenario 
is unlikely to result in low metallicity in the centers of 
these galaxies, as these dynamical effects are expected 
to have a stronger impact on the gas that is external to 
the galaxy, which is relatively more metal-poor, than the 
metal-rich gas that is located closer to the center. 

Another contrasting possibility is that the interaction 
with the nearby galaxy induces gravitational torques that 
drive metal-deficient gas from the outskirts of galaxies 
into the centers (e.g., Milros & Hernquist 1996). Sup- 
port for this idea has been seen in local interacting galax- 
ies where the radial metallicity gradient is unusually flat 
(Kewley et al. 2006; Rupke et al. 2010). This would ex- 
plain (1) the deficiency of metals, (2) the excess of nearby 
galaxies, and (3) the high SFRs and SFR surface densi- 
ties (see Section 4.6) in our galaxies. 

Our data do not decisively distinguish between these 
models without dynamical information. In particular, 
IFU spectroscopy of the ionized gas would reveal the 
presence of inflowing or outflowing gas. 

4.6. Star Formation Surface Densities 

In addition to subtracting of the continuum to identify 
nearby companions, the NB/IA excess flux images allow 
us to determine their SFR surface density, Esfr- These 
images are first converted to an emission-line luminos- 
ity, then integrated across the galaxy to determine a to- 
tal Ha or H/3 observed luminosity, and finally corrected 
for dust attenuation (Section 4.1). Six of our galaxies 
(MMT#04, #06, #07, #09-#ll) have NB imaging or 
spectroscopy that observes Ha. In addition, three other 
galaxies (MMT#13, #14, and Keck#4) have H/3 mea- 
surements. The remaining galaxies (MMT#01, #02, 
#08, #12, Keck#l, #2, and #5) are [O m]-selected. 
To obtain an estimate of the H/3 luminosity for the lat- 
ter galaxies, we use the integrated [O ill] /H/3 ratio from 
our spectroscopy (see Table 3). Finally to obtain the Ha 
luminosities, we scale the de-reddened H/3 measurements 
by a factor of 2.86. We integrate the source flux over a 
region that is 7r(FWHM/2) 2 , which is converted to phys- 
ical distance using the angular diameter distance (1" = 
5.4 kpc at z = 0.4 and 7.8 kpc at z = 0.85). This mea- 
sured area is larger due to the degradation from seeing. 
To correct for it, we determine the effective FWHM in 
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quadrature: FWHM e 2 8 = FWHM 2 - FWHM 2 eeing . Due 
to the slightly extended nature of our galaxies, this cor- 
rection to the area is typically a factor of 2 or less. We 
find a wide dispersion in the E§fr between 0.005 and 5 
M 0 yr -1 kpc -2 with an average (median) of 0.5 (0.04) 
M 0 yr -1 kpc -2 . Compared to local galaxies (Kennicutt 
1998b), our Esfr’s are on average an order of magni- 
tude higher than normal spirals, and at the low end of 
IR-selected circumnuclear starbursts. 

5. DISCUSSION 

5.1. Our [O III] X ]363 Sample 

By construction, the 20 galaxies presented in this pa- 
per were selected to be extreme. First, our SDF follow-up 
spectroscopy concentrated on ELGs identified from our 
NB/IA imaging. Then in this paper we have selected 
only a few percent of these spectra, where the lines were 
strong enough to include detectable [O ill] A4363. Half of 
our galaxies were solely targeted because of the strength 
of their [O ill] AA4959,5007 doublet (see Table 1). This 
favors emission from more highly ionized gas. Although 
the [O ill] A4363 line is occasionally seen in LINERs, we 
only found one of them based on various emission-line ra- 
tios (see Section 3.1). The remaining 19 are all extreme 
starburst galaxies, which tend to be overlooked by tra- 
ditional broad-band photometric selections. Aside from 
deep narrow-band imaging, the only other windows on 
galaxy evolution which provides access to many extreme 
SFR galaxies are extensive spectroscopy with as little 
pre-selection as possible, preferably none (such as slitless 
spectroscopy; Atek et al. 2010) or selecting galaxies with 
unusual colors (Cardamone et al. 2009; van der Wei et 
al. 2011). 

Thus these 20 galaxies lie far above the established 
correlation of star formation rate with stellar mass. In 
other words, their sSFR is between two and three orders 
of magnitude larger than the average observed in the 
local universe, or even than what is “normal” for galaxies 
at z = 0.4-0. 9. These [O in] A4363-detected galaxies 
are extremely rare in the local universe, making up only 
0.05% of those with SDSS spectra. 

Since sSFR is the inverse of the star production time- 
scale, continuing at their high observed SFRs, our ex- 
treme galaxies would produce their entire stellar contents 
in only 1%-10% of the time it takes in normal galaxies. 
This strongly suggests that we are observing them in a 
highly atypical evolutionary phase, which could last 100 
Myr or even less. If this is a phase that most galaxies 
pass through occasionally, then the ‘duty cycle’ of this ex- 
treme burst of star formation would be only ~1%-10%. 

5.2. Selection Effects 

The NB technique has two observable limitations: (1) a 
minimum EW excess at the bright end (<24 mag), and 
(2) an emission-line flux limit at the faint end. In the 
former case, our selections are limited to galaxies with 
rest-frame EWs of EWo([0 ill]) = llA 17A (z = 0.4- 
0.85) and EW 0 (Ha) ~ 20A (z = 0.24-0.40). With a 2.5cr 
emission-line flux limit of approximately 4.5 x 10 -18 erg 
s -1 cm -2 , the luminosity limit corresponds to 3 x 10 39 
and 2 x 10 4 ° erg s -1 at z = 0.40 and 2 = 0.84, respec- 
tively. These EW and luminosity limits are illustrated in 
Figure 3, and are compared to the ELGs in the SDSS. 


It can be seen that our NB imaging cannot detect the 
[O ill] emission from many SDSS galaxies, simply be- 
cause of the lower equivalent widths. For example, only 
16% of the full SDSS galaxy sample would be identi- 
fied as NB704 excess emitters at z = 0.40. However, 
this fraction increases to 90% for ELGs with moderate 
line strengths, EWo([0 ill]) > 11 A. For Ha, the major- 
ity (62%) of the SDSS would be detected at z = 0.4 with 
our NB921 imaging, and almost all (96%) of galaxies with 
EWo(Ha) > 22 A. Thus, our NB survey is biased against 
weak emission lines. We note that the lower success of 
detecting [O in] in SDSS ELGs at z = 0.4 compared to 
Ha is due to the higher metallicities in local galaxies. In 
metal-rich systems (i.e., lower T e ), 0 + is a more effective 
coolant for the ISM than 0 ++ , resulting in flux ratios of 
[O n]/Ha « 0.5 and [O m]/Ha « 0.1. Thus, our [O ill] 
NB selection is less sensitive to metal-rich galaxies. 

In addition to our general NB selection, the primary 
restriction for our sample is an [O ill] A4363 detection. 
First, this limits the sample toward high-EW galaxies, 
which generally have high sSFRs. Second, the detection 
alone will further bias our sample against moderately 
metal-poor galaxies (12 + log(0/H) ss 8.0). For exam- 
ple, in our spectroscopic survey, we have identified an 
additional ~20 galaxies that possess strong [O ill] emis- 
sion (EWo > 50A). These galaxies, however, have very 
weak detections of [O ill] A4363 (1 -3 ct). 

While an [O III] A4363 detection restricts our ELG sam- 
ple, we find that our 20 galaxies have many similarities 
to our more general ELG sample at z ~ 0.4-1. For ex- 
ample, the median (average) stellar mass of our spectro- 
scopic sample of over 200 galaxies is log (M*/M 0 ) = 8.9 
(8.7). These values are moderately higher (0.4 dex) than 
our [O Ili]-A sample. Also, as illustrated in Figure 4, our 
[O iii]-A sample shows similar [O iii]/H/? ratios, which 
is a measure of the ionization. This demonstrates that 
ELGs selected from NB surveys at z = 0.4-1 all have 
higher ionization by rs 0.2-0. 5 dex compared to local 
SDSS galaxies with similar stellar masses. 

5.3. Similarities and Differences to Green Peas and 
Luminous Compact Galaxies 

The SDSS has found rare populations of strongly star- 
forming galaxies called “Green Peas” (Cardamone et al. 
2009) and “Luminous Compact Galaxies” (LCGs) (Izo- 
tov et al. 2011). The former population was identified 
by unusual optical colors due to significant contribution 
of very strong nebular emission lines ([O ill] and H/3) 
in the r-band (hence they appear “green”). Ancillary 
data indicated that these luminous (r < 20.5) galaxies 
are at z ~ 0.1-0.35, have stellar masses of ~10 8 ' 5 -10 10 
M 0 and SFRs of ~10 M 0 yr -1 , are relatively metal-rich 
12 + log(0/H) ~ 8.7 22 , and are very rare (2 deg -2 ). Be- 
cause of the unusual color selection, selection effects are 
significant for this sample. 

Analogously, LCGs were identified from the SDSS DR7 
spectroscopic sample to have EW(H/3) > 50A, at least 2a 
detections of [O ill] A4363, and appear compact or unre- 
solved in the SDSS images. These galaxies have redshifts 
of z = 0.02-0.63 (most of them are at z < 0.35), median 
stellar masses of ~10 9 Mq, median SFRs of ~4 Mq yr -1 , 

22 Izotov et al. (2011) find that these strong-line metallicities are 
overestimated by «0.5 dex. 
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and oxygen abundance of 12 + log(0/H) ss 8.2. Since it 
is believed that green peas are a subset of LCGs (Izotov 
et al. 2011), we will only compare our sample against the 
LCG population. 

LCGs show many similarities and differences to our 
galaxies with [O ill] A4363 detections. First, both sam- 
ples occupy a distinct region in BPT and MEx diagnos- 
tics diagrams with large [O iii]/H/3 and low [N il]/H a 
(or stellar mass). While LCGs have low stellar masses 
and high sSFRs, they are on average ^0.5 dex more mas- 
sive than our sample and with lower sSFR by ^0.5 dex. 
These LCGs can be found with stellar masses above 10 10 
M 0 , while our most massive galaxy has M* « 3 x 10 9 
Mq. In addition, these galaxies are relatively metal-rich 
by ~0.4-0.5 dex compared to our sample. Also, 90% 
of our [O m]-A sample is found above z = 0.4, while 
the majority of LCGs are identified at z < 0.35. There- 
fore, while LCGs are low-mass, metal-poor galaxies, our 
[O iii]-A sample appears to be an extension of LCGs 
toward lower masses, lower metallicity, and higher red- 
shift. These differences might suggest that metal-poor 
ELGs identified with deep NB imaging may evolve into 
LCGs. Finally, in terms of surface density, it is estimated 
that there are ~20 LCGs per deg 2 . The SDF has already 
detected [O hi] A4363 in 20 galaxies over 0.25 deg 2 , and 
more of them are expected to be found (see Section 5.6). 
The higher surface density is because our survey extends 
toward lower mass galaxies, which are more common, 
and because higher sSFR are seen for galaxies at higher 
redshift (see Section 5.6). 

5.4. Comparison with Other [O ill] X 4363 Studies 

While these metal-poor galaxies are rare, there have 
been numerous efforts to identify them. Here, we com- 
pare our sample to previous [O ill] A4363 studies. Two 
of the first studies that probed a large number of low- 
mass galaxies in the local universe were Lee et al. (2004) 
and Lee et al. (2006). First, Lee et al. (2004) identified 
and detected [O ill] A4363 for 24 local galaxies from the 
KPNO International Spectroscopic Survey. These galax- 
ies, spanning Mb = -15 and -19, were pre-selected to 
have gas-phase metallicities (determined from strong-line 
methods) to be below 12 + log(0/H) = 8.2. The sample 
provided the first measurement of the MZR using the 
“direct” method. Targeting 25 dwarf irregular galaxies, 
Lee et al. (2006) extended the MZR toward lower stel- 
lar masses. Since then, more extensive spectroscopy has 
been conducted targeting low-luminosity galaxies within 
11 Mpc (Berg et al. 2012), and mining the SDSS (Brown 
et al. 2008; Izotov et al. 2006a, 2012a). 

In addition, a few other studies have succeeded at ex- 
tending the search toward higher redshift where the red- 
shift evolution of the MZR implies that more metal-poor 
galaxies should be found in the early universe. Hoyos et 
al. (2005) was the first study to detect [O III] A4363 in 17 
galaxies at z = 0.51-0.86 with gas metallicity rang- 
ing from 12 + log(0/H) = 7. 8-8. 3. These galaxies were 
selected from the DEEP2 Survey (Davis et al. 2003) and 
TKRS (Wirth et al. 2004), which are luminosity-limited 
to I = 24.5 and R = 24.3 mag, respectively. Then 
Kakazu et al. (2007) conducted follow-up spectroscopy 
of ultra-strong ELGs selected from NB imaging. While 
some level of detection for [O ill] A4363 was available in 


17 galaxies, only six had >3a [O ill] A4363 detection 23 . 
The rest were too weak, with a median detection of 

2cr. Among those with robust detections (>3a), their 
most metal-poor galaxy has 12 + log(0/H) = 7.36^0 is 
(la). They have since extended their [O ill] A4363 sam- 
ple to a total of 31 galaxies at z = 0.40-0.85 (Hu et al. 
2009). About three-fourths (23/31) of their sample are 
detected above 3 a with 9 galaxies above 5 a. Their most 
metal-poor galaxies have 12 + log(0/H) = 6.97±0.17 and 
12 + log(0/H) = 7.25T0.03 (lex). 

In addition, Atek et al. (2011) identified a sample of 
high-EW ELGs at z = 0.35-2.3 from space-based grism 
spectroscopy. With higher spectral resolution follow-up 
observations, they detected [O ill] A4363 to determine a 
gas-phase metallicity of 12 + log(0/H) = 7.5 ± 0.1 for a 
galaxy at z = 0.7. 

For a complete comparison, we compile available infor- 
mation for these studies. For the Berg et al. (2012) sam- 
ple, dust-corrected Ha SFRs are determined from the Ha 
luminosities reported in Kennicutt et al. (2008), adopt- 
ing the extinction determinations of Lee et al. (2009), and 
assuming a Chabrier IMF. We also correct their reported 
stellar masses from a Salpeter to Chabrier IMF. For the 
Hoyos et al. (2005) sample, we limit their sample to the 
seven galaxies where the [O il] line is observed. Stellar 
masses for their sample were obtained from Bundy et al. 
(2006), where a simple cross- matching of sources yielded 
only one match. SFRs reported by Brown et al. (2008) 
and Lee et al. (2004) are also corrected from a Salpeter 
to Chabrier IMF. For the samples of Kakazu et al. (2007) 
and Hu et al. (2009), we are only able to report the red- 
shifts and gas-phase metallicity for the former. In their 
final sample, only metal abundances and luminosities are 
immediately available. 

Since the rest-frame R-band absolute magnitude, Mb, 
is commonly reported, we compute it for our galaxies 
using our emission- line corrected SED. The Mg magni- 
tudes are determined as: 

M b = B + 2.5 log (1 + z) + 5.0 log (d L /10pc), (14) 

where B is the apparent magnitude at 4450A determined 
from interpolating between adjacent broad-band filters. 
These magnitudes are reported in Table 6. 

In Figure 10, we illustrate the redshifts, SFRs, stel- 
lar masses, specific SFRs, and Mg’ s for [O ill] A4363- 
detected galaxies. Our sample surveys a redshift domain 
that is similar to the studies of Hoyos et al. (2005) and 
Hu et al. (2009). In addition, our galaxies span Mg 
between -17 and -21 (excluding MMT#10, our lowest 
redshift galaxy), and have SFRs that are similar to local 
dwarf irregulars and 1 L* galaxies. 

Comparing our metallicities against their masses or lu- 
minosities, we find that half of our sample follows the 
Mb~Z or M+-Z relation of local and higher redshift 
galaxies. The remaining galaxies are located system- 
atically below these relations by 0.2-0. 3 dex. These 
galaxies have lower significance detections (3. 1-3. 6a) of 
[O ill] A4363. So while they deviate from these relations, 
the are consistent within the 95% measurement uncer- 
tainties. Deeper follow-up spectroscopy is needed to de- 
termine if these galaxies are extremely metal-poor or just 

23 One of these six does not have metallicity determination, since 
[O II] was not observed. 
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moderately metal-poor. 

Our three most metal-deficient galaxies are MMT#13 
0 = 0.47), Keck#2 {z = 0.62) and Keck#6 (z = 0.82) 
with 12 + log(0/H) = 7.53t“;“ 7.54±°;", and 7.24±° 0 f 0 , 
respectively. While these galaxies have similar metallici- 
ties (within the errors) to the most metal-deficient galax- 
ies in the local universe, I Zw 18 and SBS0335 (see Fig- 
ure 10), there are some notable differences seen in their 
emission-line fluxes. 

First, I Zw 18 and SBS0335 have [O iii]/H/3 
([O m]/[0 Ii]) ratios of 2.2 and 3.2 (9.8 and 15.5), re- 
spectively. Keck^6, which has a more robust detection 
(4.3cr) of [O ill] A4363, has a higher [O iii]/H (3 ratio by a 
factor of 1.3-2, but similar to what is seen for green peas 
and LCGs (see Section 5.3). Second, the [O ii] line is 
not detected, yielding an observed [O m]/[0 ii] ratio of 
>80. This non-detection is due to the high ionization and 
low metallicity of the ISM. In these circumstances, the 
excitation energies shift nearly completely toward 0"*" + , 
resulting in extremely low [O ii] emission. For Kecked, 
estimates using the strong-line R 23 and O32 flux ratios 
also suggest that the metallicity of this galaxy is very 
low: 12 + log(0/H) = 7.65 using the Kobulnicky & Kew- 
ley (2004) calibration 24 with a high ionization parameter, 
log(g/cm s" 1 ) > 8.8 (log (U) > -1.7). 

These galaxies are rare, but other studies have found 
them. For example, Kakazu et al. (2007) identified a 
few galaxies where [O Ii] is undetectable but very high 
S/N detections of [O ill] exist. There are also cases 
where both lines are detected, yielding [O m]/[0 11] flux 
ratios of «10-60. In addition, rest-frame optical spec- 
troscopy of Lyct emitters at z = 2.2 have measured high 
[O ill] / [O 11] ratios (Nakajima et al. 2013), which yield 
high ionization parameter estimates. 

Furthermore, several studies have shown that high-2: 
star-forming galaxies are offset from local star-forming 
galaxies in the BPT diagram with systematically higher 
[O m]/H/3 ratios (e.g., Hainline et al. 2009; Rigby et al. 
2011). A recent comparison between theoretical models 
and emission-line measurements of z = 0.5-2. 6 galaxies 
by Kewley et al. (2013) suggests that the ISM conditions 
at higher redshifts are far more extreme with high ioniza- 
tion and density. These conditions are more analogous to 
those seen in dense, clumpy H II regions of local starburst 
galaxies. As illustrated in Figure 4, our [O iii]-A sample 
is offset by 0.2-0. 5 dex in [O iii]/H/? from the average of 
local galaxies with the same stellar mass. In addition we 
find high electron densities of n e = 70-600 cm -3 . These 
observables suggest that XMPGs at z = 0.4-1 are similar 
to typical z > 1 galaxies, and can be used as a powerful 
tool for studying the ISM conditions in the early stages 
of galaxy formation. 


in the metallicity at high masses and a steep decline at 
lower masses, can be explained by a simple model where 
massive star formation enriches the ISM through stellar 
feedback. However, such enrichment can drive metals 
beyond the ISM in low-mass galaxies, and thus lowering 
their effective metallicity (e.g., Dave et al. 2011). 

Studies in the past have been mostly limited toward 
more massive galaxies (>10 9 Mq), and therefore whether 
the MZR extends toward lower masses or not is unclear 
with existing limited samples. In Figure 11(a), we illus- 
trate the MZR of AM13, which stacked SDSS spectra 
to detect [O ill] A4363 and determine metallicity using 
the same method for easier comparison. Here, the MZR 
is fitted with a logarithmic functional form devised by 
Moustakas et al. (2011): 


12 + log(0/H) 


12 + log(0/H) asm — log 



Mto 

~m7 


7 "I 


(15) 

where 12 + log(0/H) asm is the asymptotic metallicity at 
high masses, Mto is the turnover mass, and 7 controls 
the slope of the MZR at low masses. 

We note that while there are other determinations 
of the MZR using different strong-line diagnostics that 
do not require spectra stacking, there are nearly 1 
dex differences in metallicity determinations using the 
same dataset (Kewley & Ellison 2008). For this rea- 
son, we chose to use measurements obtained with the 
same method. Since AM13 used a different T e ([ O 11] )— 
T e ([0 III]) relation to determine 0 + /H + , we adopt the 
same relation, = 0.7^3 + 0.1, 25 for direct comparisons. 
We note that this difference only raises metallicity on 
average by 0.05 dex, with «0.2 dex as the largest correc- 
tion. 

We find that 9 of our 19 galaxies follow the local T e - 
based MZR to within ±2er. However, it appears that 
our entire sample has systematically lower metallicities 
at all masses. This may be the result of selection ef- 
fects, since more metal-rich galaxies are less likely to have 
[O in] A4363 detections. 

Recently, some have argued that the MZR is in fact a 
projection of a more “fundamental” plane between stel- 
lar mass, metallicity, and SFR (i.e. , M+-Z - SFR relation; 
Lara-Lopez et al. 2010). One such projection, called the 
“Fundamental Metallicity Relation” (FMR; Mannucci et 
al. 2010), suggests that galaxies with higher sSFR have 
lower metallicities compared to those at similar stellar 
masses. Efforts to describe the FMR have defined a plane 
in which a tighter correlation exists by considering a com- 
bination of stellar mass and SFR: 


AG = log 



a log 


SFR \ 

AW 1 ;’ 


(16) 


5.5. Dependence on Metallicity with Galaxy Properties 

It has been well established that the metal abun- 
dance of galaxies correlates with their stellar mass (e.g., 
Tremonti et al. 2004). The presence of such a correla- 
tion and the lack of significant scatter have provided a 
key constraint in modeling the evolution of galaxies. For 
example, the shape of the MZR, which shows a turnover 

24 This calibration yields higher metal abundances by more 0.2- 
0.3 dex against other strong-line calibrations (Kewley & Ellison 
2008). 


where a is the coefficient that minimizes the scatter 
against metallicity. We illustrate the AM13’s determi- 
nation ( a = 0.66) in Figure 11(b). A comparison against 
our sample reveals again that roughly half of our galax- 
ies follow the FMR; however, the rest of our metallicities 
fall below it. Even more interesting, some of our galaxies 
(i.e., Keck#l, #3, and MMT#=14), which followed the 
MZR are no longer within +2cr of the T e -based FMR. 
In these galaxies, the [O III] A4363 detections are robust 

25 For *3 > 2.0, we set t2 = 1.5. 
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Fig. 10. — A comparison of the properties of our [O ill]- A sample (circles) to published direct metallicity studies. The top and middle 
rows show the gas-phase metallicity vs. (a) redshift; (b) specific SFR; (c) stellar mass; and (d) rest-frame B - band absolute magnitude, 
Mb- The bottom row shows the SFR vs. (e) stellar mass, and (f) Mb- Our sample is shown as black circles while the samples of Lee et 
al. (2004), Hoyos et al. (2005), Brown et al. (2008), Kakazu et al. (2007) and Hu et al. (2009), Atek et al. (2011), Berg et al. (2012) and 
Izotov et al. (2012a) are shown as red triangles, red squares, blue diamonds, purple triangles, blue squares, red crosses, and blue crosses, 
respectively. Our only LINER (MMT#03) is denoted by orange crosses. Measurements for I Zwl8 and SBS0335 are also overlaid for 
comparison. We provide 95% C.L. for our metallicity determinations, as well as for published results. Uncertainties on stellar masses, SFRs 
and other properties are 68% C.L. Where applicable, we apply corrections to SFRs and stellar masses to a consistent Chabrier IMF. 


(S/N = 11.1, 4.6, and 7.6, respectively), and their SFRs 
are high. If anything, we see a positive or flat correlation 
between metallicity and sSFR (see Figure 10(b)), con- 
trary to what is expected from the FMR. This might sug- 
gest that a M*-Z-SFR. relation as seen for local galaxies 
may not hold in low-mass galaxies at z = 0.4-0. 9, and 
that there is a larger scatter compared to z ~ 0. 

5.6. Number Densities of 12 + log(0/H) < 7.65 
Galaxies 

With a statistically meaningful sample of [O ill] A4363 
detections, we can infer the space density of XMPGs 
from the largest emission-line sample at intermediate 
redshifts. As illustrated in Figure 3, we find that our 


[O ill] A4363 sample has rest- frame EWs of EWo(Ha) > 
100A and EWo([0 ill]) > 50A. Limiting galaxies to these 
EWs, we find a number density at z ~ 0.45 of 4.5 x 10~ 3 
Mpc -3 for Ha and 5.1 x 10 -3 Mpc -3 for [O III]. Our 
spectroscopic follow-up sampled a broadly representa- 
tive range, reaching down to EWo = 20A, and targeted 
about ~25% of the high-EW galaxies. Among these 
spectra, 16% and 11% of the Ha-selected and [O De- 
selected galaxies yielded [O III] A4363 detections, respec- 
tively. Extending these statistics to the full high-EW 
sample, we estimate that the expected number density 
of galaxies with [O in] A4363 detections at z = 0.3-0.65 
is 7 x 10 -4 Mpc -3 for Ha-selected galaxies and 6 x 10 -4 
Mpc -3 for [O m]-selected galaxies. This corresponds to 
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Fig. 11. — Comparison of T e -based gas metallicity against (a) stellar mass, and (b) a projection of the “FMR” with fi = 0.66 for our 
[O III]- A sample (black circles). The MZR and FMR of the SDSS derived using the T e method by AM13 are shown in grey with dashed 
lines indicating zb2<r. Here, the scatter of AM13 is not the true intrinsic scatter from individual galaxies. Rather, it reflects the dispersion 
for stacked spectra in various M*-SFR bins. This cr is likely to be larger than the intrinsic scatter because it is weighted more toward 
high-SFR outliers (fewer galaxies are available in these bins; B. Andrews 2013, priv. comm.). Metallicity uncertainties are reported at 95% 
C.L., while 68% C.L. are reported for mass and SFR uncertainties. For direct comparisons, we adopt the same £ 2^3 relation as AM13. 


a surface density of one per 13 arcmin 2 (H a; z = 0.3- 
0.5) and one per 8 arcmin 2 ([O III]; z = 0.3-0.65). About 
one- fifth of these galaxies will have 12 + log(0/H) < 7.65. 

The higher rate of XMPGs («1%) among star-forming 
galaxies at high redshift compared to the local universe 
(recall 1 in 1700 SDSS galaxies are XMPGs) could be a 
reflection of galaxy evolution. In particular, several lines 
of evidence have shown that the SFRs and the sSFR’s 
are an order magnitude higher at z > 1 than the local 
universe (e.g., Noeske et al. 2007; Ly et al. 2011b, 2012b; 
Lee et al. 2012). Future surveys at z > 1 with JWST 
and “WFIRST” should find many more XMPGs sim- 
ply because: (1) the M+-Z relation evolves toward lower 
metallicity at high redshift (see e.g., Erb et al. 2006); and 
(2) the epoch of where the SFR density peaked (Hopkins 
2004; L07), and where most of the mass assembly oc- 
curred (Dickinson et al. 2003; Rudnick et al. 2003), is at 
2 ~ 1-3. 

Morales-Luis et al. (2011) reported that the number 
density of XMPGs in the local universe from the SDSS 
is ~ 1.3 x 10 -4 Mpc -3 , which is similar to what we find in 
our study. This comparison would suggest that no evolu- 
tion is seen in the number density of XMPGs. However, 
we caution against making such a comparison since the 
local estimate is highly uncertain due to a statistically 
small sample that suffers from incompleteness and com- 
plex selection functions. 

6. CONCLUSIONS 

We have conducted an extensive optical spectro- 
scopic survey of nearly 900 ELGs identified in the 
SDF from NB or IA imaging with Keck/DEIMOS and 
MMT/Hectospec. Upon inspecting our spectra, we 
found that we have 20 galaxies where the weak auro- 
ral [O in] A4363 line is detected at >3 ct. These detec- 
tions allow us to determine the electron temperature of 
the ionized gas, and directly determine gas-phase metal- 
licity, using the T e method. Excluding our one object 
which we believe is a LINER, we find that the remaining 


19 galaxies are metal-poor with a median 12 + log(0/H) 
of 7.75 (0.13 Zq), and the most metal-rich galaxy having 
12 + log(0/H) = 8.40 (0.5 Zq). We find that 4 of our 
galaxies are XMPG (12 + log(0/H) < 7.65). 

Our most metal-deficient galaxy with a >4 a detection 
of [O in] A4363 has an oxygen abundance (12 + log(0/H) 
= 7.241 q 3 q) that is similar to I Zw 18 and SBS0335. This 
galaxy lacks an [O n] detection, which is caused by a high 
ionization parameter and low metallicity that drive most 
of the oxygen species to be in a doubly ionized state. 
In addition, we find that all of our galaxies have high 
ionization, as measured by [O iii]/H/3, and high electron 
density (n e = 70-600 cm -3 ). While this is not typi- 
cal for SDSS ELGs, there are rare populations of highly 
star-forming low-z galaxies found in the SDSS that have 
high ionization. In addition, these properties are seen for 
strongly star-forming z > 1 galaxies. This suggests that 
XMPGs at z ~ 0.4- 1 are analogs to high-z star-forming 
galaxies in terms of their ISM properties. 

Our exceptionally deep and complete multi- wavelength 
dataset in the SDF also allows us to determine other key 
properties of these galaxies. What we found is that these 
galaxies are exceptional in several ways. 

First, with stellar masses between 5 x 10 6 and 3 x 10 9 
M q (with a median and average of 3 x 10 8 and 2 x 10 8 
Mq, respectively), we find that our galaxies have high 
specific SFRs, which suggests that these galaxies have as- 
sembled their stellar contents in only ~100 Myr. Second, 
their vigorous star formation is highly concentrated into 
a small central starburst of high surface density (an aver- 
age of Esfr = 0.5 Mq yr -1 kpc -2 ). Third, they have an 
excess probability of having nearby companions (within 
a projected 100 kpc), which we have determined to have 
strong emission lines at a similar redshift. The presence 
of nearby companions in two-thirds of our galaxies sug- 
gests that galaxy-galaxy interactions may be responsible 
for the metal deficiency and the observed high sSFRs 
and high SFR surface densities. In particular, outflow- 
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ing or inflowing gas can occur through stellar feedback 
that drives metal-rich gas out of the ISM or gravitational 
torques that can cause metal-poor gas to lose angular 
momentum and fall into the center of the galaxies. 

Some researchers have proposed that the M+-Z re- 
lation must be extended to include a third parameter, 
the SFR (Lara-Lopez et al. 2010; Mannucci et al. 2010). 
Since these relations predict lower metallicities for galax- 
ies with unusually high sSFRs, about half of our galaxies 
are indeed consistent with the FMR. However, there are 
two inconsistencies that we find. First, our most star- 
forming galaxies (3 of them) have higher metallicity than 
what is expected if they followed the FMR. This sug- 
gests a flatter or positive correlation between sSFR and 
Z , contrary to an inverse correlation found for the FMR. 
Second, it appears that our sample has systematically 
lower metallicity by ?»0.2 dex than the FMR. This re- 
sult is only robust at 2cr confidence, and we suspect that 
selection effects (i.e. , required detection of [O ill] A4363) 
may be biasing our sample toward lower metallicity. 

The failure of any simple relation to predict O/H can 
be seen by the large scatter of our metallicities versus 
stellar mass or sSFR in Figure 11. The large intrin- 
sic scatter overwhelms any other trends with metallicity 
which might be present. For example, our sample con- 
tains 19 galaxies with redshifts ranging from z=0. 38-0.88 
(excluding the lowest redshift galaxy), a range of look- 
back times from 4 to 7 Gyr. This period of cosmic his- 
tory is known to span a great deal of evolution in many 
properties of the overall galaxy population as a whole. 
Nonetheless our galaxy sample shows no trend of metal- 
licity with redshift. Across all of these redshifts, we find 
galaxies with half solar abundances mixed in with very 
metal-poor galaxies. We believe that this large cosmic 
scatter of abundances is larger than our measurement 
uncertainties, and is genuine. 

The DEIMOS data presented herein were obtained at 
the W.M. Keck Observatory, which is operated as a scien- 
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TABLE 6 

Luminosities, Stellar Properties, and Star Formation Rates 
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Note. — All values assume a Chabrier IMF, and where applicable, dust attenuation follows COO formalism. 68% confidence measurement 
uncertainties are reported. 

a One-fifth solar metallicity results are given Cols. (7), (9), and (11) while those obtained with solar metallicity models are provided in Cols. 
(8), (10), and (12). 

^This source is likely a LINER, so SFR estimates are not trustworthy. 

c Since the Ha line falls at the edges of the filters, these Ha SFRs are affected by accurate tracing of the NB filter profiles. 
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APPENDIX 

A. ACCURATE FLUX CALIBRATION OF MMT/HECTOSPEC AND KECK/DEIMOS SPECTRA 

Since the spectroscopic data were obtained from two different instruments, MMT/Hectospec and Keck/DEIMOS, 
using different observational configurations (fibers versus slits) and in various observational conditions (seeing: OU 5— 
1"5), the flux calibration of our data is problematic. To better unify the spectroscopic data, we use the following 
approach. 

First, spectro-photometric standard stars were observed. These observations were generally obtained on the same 
night; however, in some cases, these calibration data were taken a few days apart. Reducing the calibration data in 
identical or similar ways to our science spectra, we determine the sensitivity function(s) using standard iraf processing 
techniques. We then apply these sensitivity functions to our data to yield “first-pass” flux-calibrated spectroscopic 
products. 

Second, we examine the reliability of our flux calibration for each slit mask or fiber configuration by comparing them 
to broad-band photometric data. Here we only consider the brightest galaxies such that the continuum is well-detected. 
For the MMT and DEIMOS data, we restrict the calibration sample to Rc < 22.0 mag and i' < 23.0 mag, respectively. 
This typically has 10 galaxies per set-up. To avoid the effects of OH skylines, we generate smoothed spectra where 
a boxcar median is used. From there we convolved our smoothed spectra against the filter responses and determine 
fluxes. These fluxes are compared against the photometric data to examine slit/fiber losses and if any second-order 
wavelength-dependent corrections are needed. In most cases, the comparison against photometric data only showed a 
systematic offset due to slit/fiber losses. This comparison showed that the rms is typically between 10%-20%. 

Next, to ensure that the flux calibration is reliable, we use our measurements of nebular emission lines from NB 
imaging data to compare against the spectroscopic measurements. The NB fluxes are derived from Equation (6) but 
differs for each NB filter. We use the spectroscopic redshift to correct the NB flux for when the emission line falls 
in the filter’s wing. This comparison is illustrated in Figure Al. For 116 NB816 emitters with MMT spectroscopy, 
we find excellent agreement with a small offset (-0.036 dex) and an rms of 0.155 dex. Similar results are also seen 
for 108 NB704 emitters. We note that these calibration results are consistent with those obtained by the Hectospec 
instrument team (Fabricant et al. 2008). In addition, the DEIMOS fluxes also show great agreement for 58 NB921 
emitters. Similar results are also seen for 9 NB704, 39 NB816, and 8 NB973 emitters. This consistent agreement at 
multiple wavelengths supports our wavelength-dependent flux calibration. Note that in our second step we only use 
broad-band data, so this comparison against NB data is an independent test. 

Finally, a second test that we conducted is to examine the emission-line flux determinations for the same galaxies in 
different observing conditions. This was only feasible with MMT because of the wide field-of-view to target multiple 
galaxies simultaneously. In total, we have 28 galaxies with two spectra. We illustrate in Figure A2 the determination 
of fluxes for various emission lines that span rest-frame wavelengths of 3700A-6730A and observed wavelengths of 
4000A-8300A. We find that between different MMT setups, the same emission-line fluxes are obtained with an rms of 
0.15-0.20 dex across a wide spectral range, which further supports our approach for accurate flux calibration. 

We note that the effects of differential refraction are mitigated for the MMT spectra, as an atmospheric dispersion 
compensator (ADC) is used and is known to eliminate issues affecting the flux calibration of the spectra (Fabricant et 
al. 2008). DEIMOS, however, does not have an ADC. The wavelength-dependent effects of differential refraction are 
seen in some of our DEIMOS observations when comparisons are made against broad-band data. Since our data are 
calibrated against broad-band photometry, as discussed above, these issues are resolved (to first order). 

B. REST-FRAME OPTICAL SPECTRA FOR [O hi] A4363 GALAXIES 

Here, we provide the rest-frame optical spectra for our sample of 20 galaxies in Figures B1-B3. These plots are 
limited to a rest-frame wavelengths of 3400A-5050A. Our figures are limited toward lower fluxes to show the weak 
emission lines. Since the strong lines are saturated in these figures, we refer readers to Table 3, which provides the 
emission-line strengths. Overlaid in blue on these spectra are the locations of OH sky-lines with the OH line strengths 
indicated by the darkness of the blue shades. We also overlay strong and weak nebular emission lines in red. We 
previously discussed in Section 2.2 the acquisition and reduction of these spectra, and will discuss the absolute flux 
calibration of our spectra in Appendix A. 
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Fig. Al. — Comparison between spectroscopic flux measurements against NB fluxes. The left (right) panel shows the MMT (DEIMOS) 
flux measurements for NB816 (NB921) emitters. Measurements for Ha+[N ii], [O ill], H/3, and [O ii] are shown as red, green, purple, and 
blue, respectively. Diamonds denote those measurements where more than one emission line enters the NB filter (e.g., Ha and [N ii], or 
[O ill] A4959,5007). These figures demonstrate that the flux calibration is precise and accurate for both MMT and DEIMOS. Other NB 
emitters have been inspected and similar results have been found. 




Fig. A2. — Comparison between emission-line fluxes for various MMT setups. The top panels show the flux comparisons while the 
bottom panels show the differences against the flux. The different colors and symbols correspond to a given emission line (e.g., purple 
squares for [O ii] AA3726,3729). Measurements that are reliable at >3cr in both spectra are marked with filled circles. The shaded grey 
regions in the bottom panels show the median and drier. The median is computed using only >3cr measurements. 
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Fig. Bl. — MMT/Hectospec spectra for seven galaxies with full spectral coverage and detection of [O ill] A4363 at >3cr. The most 
common nebular emission lines between 3400A and 5010A are denoted by the red dashed lines. OH sky-lines are indicated by blue vertical 
bands where the strength of the sky-lines is denoted by their shade of blue (darker is stronger). Grey shaded regions are those affected by 
atmospheric absorption (A- and B-bands). 
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Fig. B2. — Same as Figure B1 but for another seven galaxies observed with MMT/Hectospec. 
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Fig. B3. — Keck/DEIMOS spectra for six galaxies with full spectra coverage and detection of [O ill] A4363 at >3 cr. The most common 
nebular emission lines between 3400A and 5010A are denoted by the red dashed lines. OH sky-lines are indicated by blue vertical bands 
where the strength of the sky-lines is denoted by their shade of blue (darker is stronger). Grey shaded regions are those affected by 
atmospheric absorption (A- and B-bands). 



